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HYDROLOGIC STUDIES ELECTRONIC COMPUTERS TVA 


Willard Snyder,! ASCE 


SYNOPSIS 
Application electronic computers TVA watershed studies pre- 
sented. comprehensive multiple regression analysis program for the 
IBM 704 computer described. The use this program developing 
and studying hydrologic prediction equations illustrated discussion 


numeric results fitting two equations data. 


INTRODUCTION 


Hydrologic studies the TVAare carried under research 
and development watersheds, which are tributary areas the Tennessee 
River. The program administered Advisory Committee Tributary 
Watersheds. The purpose the program and the role hydrologic investi- 
gations are the first two statements Understanding and Beliefs, 
formulated the Advisory Committee. These are follows: 


The welfare the people the clear goal the tributary watershed 
program, and their full participation essential its success. 

watershed activities are broadly distinguishable. Both may carried 

the same watershed area, but each designed achieve different but related 

ends. Research activities are designed produce reliable facts from experi- 

ments with single defined combinations controlled practices. These facts 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper 
ofthe copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. February, 1960. 


Head, Hydrology Sect., Hydr. Data Branch, TVA, Knoxville, Tenn. 
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are basic the guidance action programs developmental watershed 
projects, where practices are not controlled. 


TVA has variety programs for development the resources the 
watersheds the Tennessee Valley. These include cooperative research 
projects with state universities, test fertilizer research, 
and reforestation activities. programs, when undertaken sepa- 
rately, not achieve maximum development the welfare 
the people the primary concern, watershed research and development must 
undertaken integrated basis. The problems for which the best possi- 
ble solutions are sought are found the interaction such factors land- 
use, soils-water relationships, and the socialand economic levels the areas. 

watershed approach integrated solution necessary one wishes 
completely evaluate program development and improvement. Starting 
with rainfall one the resources area, the changing hydrology must 
considered relation the total development resources that area. 
Therefore, hydrologic studies and analyses are undertaken one part 
major program. Causes and effectsof many varying factors must evaluated 
precisely possible watersheds designated for research activities. 
These results form basis for development programs other projects. 

Evaluation hydrologic data integrated program requires development 
ways and means for (1) the adjustment data eliminate variability in- 
herent uncontrolled experiments, (2) detection hydrologic changes, and 
(3) the synthesis hydrologic response under changing watershed practices. 
accomplish these itis necessary formulate and evaluate effects vary- 
ing amounts and patterns rainfall, soil moisture levels and rates loss, 
among others, and eliminate their effect volume runoff, peak rates 
discharge, and rates sediment movement the residual data, after this 
hydrologic adjustment, are less variable than before adjustment. These re- 
sidual data then form the basis for synthesis external effects, 
such watershed management, vegetal cover, and possibly soils types. 

Much the analytic effort involved hydrologic studies is, course, re- 
ducible established computational forms. Where numeric procedures exist 
canbe devised they can usually programmed for electronic computation. 
From practical standpoint, however, only those procedures should pro- 
grammed which the volume repetitive work justifies the initial cost 
program preparation. 


TVA COMPUTERS 


Electronic computation hydrologic studies TVA began September, 
1955 with use facilities the National Laboratory Oak Ridge, Tennessee. 
December, 1957 the TVA Division Water Control Planning obtained 
small Royal-McBee LGP-30 computer. June, 1958 installation high 
speed IBM 704 computer was completed. Hydrologic computations have been 
programmed for both TVA machines. 

Programming for computers TVA user basis. That is, programs 
for hydrologic computation are written engineers engaged hydrologic 
studies. Programs dealing with system operations are written engineers 
engaged river control activities. this end, selected engineers and engi- 


neering aides from the Hydraulic Data Branch, and every other using Branch, 
have been trained computer programming. 
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Hydrologic problems programmed for the small LGP-30 computer were 
those which arose during basic information, such computation 
discharge from time-gage height abstracts, computation daily 
suspended sediment loads from discharge and concentration data, and compu- 
tation indices antecedent precipitation. 

Hydrologic problems programmed for the large IBM 794 computer were 
those which contain complex voluminous computations. Two programs have 
been very useful. One computes daily values soil moisture index from 
rainfall, runoff, and evaporation for extended periods record. The second 


program for comprehensive multiple regression analysis based the 
method least squares. 


THE MULTIPLE REGRESSION PROGRAM 


The multiple regression program very definitely the most valuable 
the current hydrologic computer programs. The method least squares 
technique for evaluating the simultaneous effect many factors has become 
well established hydrology. Coupling this technique with regression analy- 
sis makes the scope analytic results almost unlimited, and appears offer 
the best approach hydrologic evaluation necessary part the TVA 
Tributary Watershed Studies. 

The multiple regression program used TVA was written include five 
major phases analysis fitting prediction equations observational data. 
These will described sequential order computation. 

Product Summation.—Input the program consist individual 
vations and associated parameters which have been computed for each obser- 
vation. Generation triangular matrix product-summations the first 
phase computation. This matrixis the setof coefficients the simultaneous 
normal equations for least-squares fitting. matrix 500 obser- 
vations can generated. The triangular matrix printed. 

Simultaneous symmetrical square matrix generated from 
the triangular matrix. This matrix inverted and the regression coefficients 
are computed. Both are printed. Following this, various statistical quantities 
resulting from fitting the equation are computed and printed. These include 
sums squares for regression and for residuals, the correlation coefficient, 
the error prediction and the error regression, the standard errors the 
regression coefficients, and the arithmetic mean each variable. 

Predicted Values and Prediction Errors.—After the statistical results 
fitting are printed the observations are into the computer. The pre- 
dicted value for each observational set then computed from the regression 
coefficients and the independent variables. The predicted value subtracted 
from the observed value give the prediction error for each observation. 
Both predicted value and prediction error are printed. 

Predicted Values and Variance Prediction.—To this point the re- 
gression program the operator has not been able exercise any options 
program computation output. The last two phases compu- 
tation are optional. 

sometimes desirable totest the statistical variability the regression 
equation aprediction device. The prediction errors computed for 
the initial data set not meet this requirement. However, standard statistical 
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procedures can found for this computation. This procedure, operation 
the inverse matrix, was written into the program. 

use this optional portion the program list synthetic observations 
must prepared advance, and punched oncards placedon magnetic tape 
along with the list observations used fitting. The only requirement set- 
ting the synthetic list set the independent variables values which 
desired test the variability prediction. The machine computes the 
predicted value and the variance the each synthetic observation 
the set. For simple X-Y correlation this portion the regression pro- 
gram would compute the confidence limits about the line 
regression. For many independent variables the model complex and the 
arithmetic becomes quite lengthy. More than 200 product-summations are re- 
quired. For evaluate the variance for single prediction 10- 
term equation, and the computer tremendous time saver this operation. 

Significance Variables.—The last phase the computer program was 
written touse testing the different independent variables for statistical sig- 
nificance. This phase the program also optional. 

true that the standard the can used 
form test significance the statistic This test, however, limit- 
one variable. When several variables incombination are tested for 
significance, the regression equation must re-evaluated with this combi- 
nation terms left out. The resulting change the sum squares absorbed 
the regression forms the basis for analysis variance, using the statistic 
“F” determine the significance the entire combinationof variables under 
consideration. 

The TVA multiple regression program allows any number re-evaluations 
the regression equation with any number any combination variables left 
out. The combinations terms dropped are independent the order which 
the re-evaluations are made. With this complete freedom choice, the ana- 
lyst specifies beforehand the various re-evaluations which are desired. For 
each reduced equation the program causes computation the reduced inverse 
matrix and regression coefficients and the other various elements deter- 
mined for the first evaluation. The predicted values, however, are not comput- 
ed. Comparison the solutions for all the reduced equations, with the initial 
total equation can then made all these elements the solution. These, 
course, are addition the conventional analysis variance. 


EQUATION FOR PEAK DISCHARGE 


Some examples the equations being evaluated the computer program 
for regression analysis may interest. is, perhaps, seeming contra- 
diction, that proper utilization computers requires better working know- 
ledge mathematics than when investigations are based purely empiri- 
cisms and short-cut assumptions. However, computer programs must follow 
precise part the goal computer utilization more precise 
analysis. 

One particularly vexing problem TVA hydrologic studies has been the 
predictionof peak discharges for some small, single-practice research water- 
sheds western North Carolina. The watersheds are operated coopera- 


McGraw-Hill Book Co., New 1952, 202. 
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tive project between TVA and North Carolina State College. Four watersheds, 
cally designed experiment. Results from four different covers the four 
watersheds will studied analysis variance following latin square 
design. Preliminary this analysis will necessary establish the re- 
lationship the important hydrologic and meteorological variables, that 
data can adjusted for varying meteorological conditions. 

these small watersheds the normal hydrologic procedure correlating 
peak discharge with peak-producing rain, maximum intensity, antecedent 
moisture, and one, would not produce workable peak equations. The runoff 
the measuring flumes responds quickly changes rainfall intensity 
that there are quite often several peaks during storm, each peak represent- 
ing burst rainfall. After several trials analysis, was decided that the 
patternor time distribution rainfall would have The follow- 
ing development leads such equation. 


TIME_IN 


RAINFALL 


FIG. 1.—DEFINITION RUNOFF PERIODS WITH RELATION 
PEAK DISCHARGE 


Any ordinate storm hydrographcan predicted froma unit hydrograph 
and incremental volumes runoff. Peak discharge, being one ordinate, can 
predicted equation which expresses the unit hydrograph principle: 


which peak discharge, represents the set unit hydrograph coef- 
ficients, and the volume runoff periods. 

The quantities Eq. (1) can considered time-distribution coefficients 
and they not have atuniform intervals.3 For the western 
watersheds the coefficients were assigned time definition shown Fig. 
this schematic diagram rainfall and the resulting hydrograph are plotted 
acommon time scale. Zero time taken the time peak, and the time 
definition the coefficients extends backward the beginning rainfall. 


“HYDROGRAPH ANALYSIS THE METHOD LEAST SQUARES,” 
Snyder, ASCE Proceedings Paper No. 793, September, 1955. 
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The volumes rainfall, can obtained from recording rain gage. 
necessary devise method predict runoff from rainfall for each the 
four time periods. Runoff excess rainfalland can computed for each 
period from rainfall minus loss: 


where loss considered negative. 
Loss was considered two parts. The first part was initial amount 
which varied with season but which was independent storm rainfall. The 


second part was progressive amount which accumulated during the storm 
rate governed the accumulation storm rainfall, but also dependent upon 
initial soil moisture. soil moisture index, used this relation, was com- 
puted modifications previously reported budgetary system which com- 
putes two-level index.4 The deep-level index varies slowly, and was used 
the expression for season. The upper-level index varies rapidly, and was 
used surface moisture index. This two-level index was calculated for 
each day separate computer program, using daily rainfall, daily runoff, 
and adjusted monthly pan evaporation. 


The equation used express the total loss terms the two component 
parts was 


which the loss for four time periods, represents the surface- 
layer soil moisture index the beginning each time period, the deep- 
layer soil moisture the beginning the storm, denotes the rainfall during 
each period, No, Mo, and are threshold values respective variables, and 
and are relationship coefficients. 

Eq. was used because the simplest possible expression inter- 
action effect two causative factors, plus the independent effect third 
factor.5 

The surface soil moisture index, may considered increasing di- 
rectly with rainfall. This approximately correct when direct runoff 
small portion the rainfall. From this derived: 


surface-layer moisture index beginning storm 

Eq. can used with Eq. define the loss for each the four time 
periods. Only rainfall during the preceding periods, initial moisture index 
the surface layer, and the moisture index the deep layers are required. Com- 
bining the resulting equations for period losses with Eq. these expressions 


for runoff can substituted into Eq. The resultant prediction equation for 
peak discharge algebraically complex its coefficients, but simple its 


“COMPUTATION EVAPORATION AND EVAPOTRANSPIRATION FROM METE- 
OROLOGICAL OBSERVATIONS,” Kohler. Paper prepared for presentation 
AMS Meeting Chicago, March 19-21, 1957. 
“SOME LINEAR MODELS FOR MULTIPLE REGRESSION ANALYSIS HYDROLO- 
GY,” Snyder. Duplimat Reproduction Tennessee Valley Authority, Knox- 
ville, Tennessee, May, 1956. 
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STUDIES COMPUTERS 


variables. collectirg coefficient terms the equation can simplified the 
following: 


Eq. the variables and which the peak, Qp, related have 
been previously defined. The coefficients, however, result from algebraic sim- 
plification and have real hydrologic interpretation. The numerical value 
the coefficients can determined the methodof least squares. The results 
obtained thus fitting the equation data can analyzed multiple re- 
gression analysis. 

Eq. has beenevaluated for several watershed cover data from the 
western North Carolina project. Initial results for this interim stage the 
project appear entirely satisfactory. Some numerical results are shown 
Table table multiple correlation coefficients and the conventional 


TABLE 1.—PEAK DISCHARGE REGRESSION ANALYSES FOR WESTERN NORTH 
CAROLINA COOPERATIVE RESEARCH 


Pasture 


Watershed 


Corn 


Wheat 


Number storms 


Multiple correlation 


coefficient 0.76 
0.014 0.256 


0.069 0.742 


Error regression 


Average peak 


for errors and averages are cubic feet per second. 


statistical parameter, the error the regression, are tabulated for compari- 
son with the average peaks which were observed under different covers. The 
multiple correlation coefficients show that, general, the equation making 
significant adjustments the data. However, this relative measure and 
not too useful. more absolute check the error regression. Table 
shows that the errors are small for pasture and over-grazed pasture, but that 
larger errors are associated with corn and wheat. 

Another way studying the equationis the magnitude peaks whichare 
predicted and the variance, probable variation, those predictions. Table 
shows predictions for storm high intensity with and without rainfall im- 
mediately preceding. The predicted peaks are reasonable for all covers. For 
both the short duration and long duration storms, the peaks are lowest for 
pasture, slightly higher for over-grazed pasture. However, predicted peaks 
are highest for corn for the short duration summer intense storm, and 
highest for wheat for the winter intense storm. This follows the known cover 
cycle poor cover under wheat during the winter-time, following fall planting 
and poor cover under corn inspring andearly summer following spring planting. 


Over- 
grazed 
0.42 0.79 0.86 0.90 0.88 
0.118 0.009 0.018 0.102 0.155 
0.374 0.124 0.166 0.660 1.018 


The errors prediction givenin Table are simply the square roots the 
respective prediction variances. The probable variations are given this 
form that the units are the same for the predicted peaks. The errors 
prediction show that, general, the differences between peaks for various 
cover are real. However, some improvements must still made. Variances 
for corn and wheat covers must reduced. Special attention may also have 
given Watershed This watershed, under certain storm conditions, 
large amounts subsurface flow, producing very different peaks than 
flow were from surface runoff. consequence, peak discharge data from 
this watershed are more variable than from the other areas, and the average 
residual errors also tend large. attempting improve the prediction 
peak discharges study will first made the prediction errors for the 
individual storms. These data are available from the machine printout sheets, 
and attempt will made associate the errors with cover watershed 
condition varies with the annual cultivation cycle. 


TABLE 2.—PREDICTED PEAK AND ERROR FOR WESTERN NORTH CAROLINA 
COOPERATIVE RESEARCH 


Wheat 


Predicted peak 


Standard error 
prediction 


Storm Parameters (Inches): 0.50, 1.00, M=1.00, N=20.0 


Predicted peak 


Standard error 
prediction 


are cubic feet per second. 


EQUATIONS FOR VOLUME RUNOFF 


equation for prediction volume storm runoff from the small water- 
sheds western North Carolina also under development. This equation 
very similar tothat used for prediction peak discharge. The similarity lies 
the definitionof rainfall quantities short periods time. There dif- 
ference, however, that the volume equation has the most precise time defi- 
nition the beginning the storm where loss rates are high and changing 
rapidly. 

Fig. diagrammatic definition the periods rainfall used the 
volume runoff equation. Fig. rainfall and the resulting hydrograph 
are plotted common time scale. However, zero time Fig. taken 


the beginning storm rainfall, and time definition extends forward the end 
rainfall. 


Over- 

grazed 
0.40 5.80 1.64 0.37 1.01 0.03 1.22 

0.05 1.14 0.68 0.04 0.08 0.55 0.51 

0.43 6.40 8.08 1.04 2.13 3.50 8.99 


The volume runoff equation, based the rainfall periods Fig. 
and soil moisture is: 


which surface plus subsurface runoff, are rainfall amounts dis- 


RAINFALL 


Surface 
plus 
Sub-surface Runoff 


FIG. 2.—DEFINITION RUNOFF PERIODS WITH RELATION 
STORM RUNOFF VOLUME 


tributed time periods, and and are moisture parameters described pre- 
viously. 

Only few analyses this equation can reported this time. Results 
these fittings are encouraging and are shown Table Multiple corre- 


TABLE 3.—RUNOFF VOLUME REGRESSION ANALYSES FOR WESTERN NORTH 
CAROLINA COOPERATIVE RESEARCH 


Watershed 
Pasture Pasture Wheat 


Number storms 114 


Multiple correlation 
coefficient 0.75 


Error regression 0.003 


Average runoff 0.029 


for errors and averages are area inches. 


lation coefficients are high and errors regression are low. Some study 
the residual errors these fittings indicates there may improvement possi- 
ble taking into account the volume rain above some critical intensity, 
perhaps for each rainfall period. 


0.94 0.90 
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COSTS COMPUTATION 


would impractical evaluate many equations the distributed rain- 

fall type desk calculator. The time required for solution these lengthy 
equations would allow only few trials. Many runs the multiple regression 

program have been possible onthe 704 computer. cost comparison was made 

after analyses. Computer costs for these runs totalled $663. This does 

not include cost writing the program but does include costs checking out 

the program onthe machine. was estimated that the same analyses using 

desk calculators would have cost about $10,000. This estimate based 

extrapolation costs for comparatively short equations. There allow- 

ance made this extrapolation for increasing number operator errors 

due the tedium lengthy calculations, nor the possibility increased 

personnel turn-over due pure boredom. 


OUTLOOK 


The electronic computer will used increasing rate hydrologic 
studies. The high speed computation makes possible the accomplishment 
routine work quickly and efficiently. The working time hydrologists and 
aides, released from much routine computation, can used develop new 
methods, new techniques, and new machine programs. The low unit cost 
analysis when computedelectronically allows rapid trialof these new methods 
over many observations and many data sets. 

Rapid strides are being made techniques data analysis. 
The great challenge, however, lies hydrologic synthesis. Hydrologists will 
not making full use computer capabilities until techniques are developed 
whereby proposed watershed programs canbe evaluated integrating the ef- 
fects changeson many small sub-areas. This has beendone for water yield, 
based infiltration classifications the sub-areas. But important that 
practical techniques integration developed for changes flood peaks, 
shape the hydrograph, and transport sediment. These are factors which 
can unquestionably improved good land management and should serve 

the bases for evaluating proposed programs. 

All phases computer applications, basic processing observations, 

sis data, and research for new techniques are necessary steps continu- 

ing program. The end result this program make more and better hydro- 

logic information available agronomists, foresters, and hydraulic engineers 
who are working improve the welfare the people. 
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SCALES VISCOUS ANALOGY MODELS 
FOR GROUND WATER STUDIES 


Jacob Bear! 


SYNOPSIS 


coefficient, two liquid flows, and capillary effects. 


INTRODUCTION 


known the Hele-Shaw model. 


Proceedings the American Society Civil Engineers 


Scales for viscous fluid models used for studying various problems 
flow through porous media are described. Consideration given 
some model techniques and limitations. The discussion covers scales 


for length (including distorted scales), time, discharge, volume, storage 


The movement fluid through porous media complicated phenomenon, 
and the mathematical equations that describe can solved only for few 
very simple cases. Difficulties are removed large extent when model 
studies replace the theoretical approach. Among the various types models 
which are used for studying ground water movement are the sand model, the 
electric analogy, the heat analogy, the membrane analogy, and the viscous 
flow analogy. This paper concerned with the viscous flow analogy, also 


The viscous flow analogy model consists two parallel plates mounted 


together with capillary interspace (0.5-3.5 mm) between them. The analogy 


based the similitude between the differential equations which describe 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 


ciety Civil Engineers, Vol. 86, No. February, 1960. 


Graduate Research Engr., Univ. Calif., Berkeley, onleave from Technion, Israel 


Inst. Tech., Haifa, Israel. 
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the field saturated flow water through porous media and those for lami- 


nar flow viscous fluid through the capillary interspace between two 
parallel plates. 


The first model this kind was developed Hele-Shaw2,3 1897- 
1899, for the study flow patterns around variously shaped bodies. 
1936, was the first suggest the application this model 
ground water flow studies. Since 1940, the model has been much improved, 
developed, and intensively used the Hydrological Laboratory the Govern- 
ment Institute for Water Supply, The Hague, Netherlands, for studying various 
problems ground water movement5,6,7,8,9; the Israel Institute Tech- 
nology, Haifa; and Water Planning for Israel Tel-Aviv. Other studies and 
applications have been reported 11,1 

The viscous flow analogy can used for almost any study two 
sional flow porous media whether steady not. has been used for seep- 
age through dams and from channels, for tidal fluctuations the water table, 
for ground water balances, for studies the movement the salt water- 
fresh water interface coastal aquifers, and other problems. 1956, this 
type analogy was further improved include the three dimensional case 
horizontal flow aquifer with fluctuating water 

The principles the viscous flow analogy are well known and described 
many articles dealing with this type model, yet complete description 
the computations scales for the various cases missing. This paper 
presented effort fill this gap. 


“Experiments the Nature Surface Resistance Water and Streamline 
Motion Under Certain Experimental Conditions,” Hele-Shaw, Trans. Institution 
Naval Architects, vol. 40, 1898, pp. 21-46. 

“Streamline Motion Viscous Film, Rept. 68th Meeting British Assn. for Ad- 
vancement Science,” Hele-Shaw, London, 1899, pp. 136-142. 
“Grundwasserstromung,” Dachler and Julius Springer, Vienna, pp. 118-120, 
1936. 

“Ein Modelproef Ter Bestudeering Van Nietstationaire Bewegingen Van Het Grond- 
water” model test for studying non-steady ground water flow); Dietz, Water 
(Holland) vol. 25, no. 23, The Hague, 1941. 

“Erwaringen Met Modelonderzock Hydrologie” (Experience gathered model 
research for hydrological investigations); Dietz, Water (Holland), vol. 28, no. 
The Hague, 1944, 17-20. 

“Modele pour des problemes écoulement l’eau souterraine employé 
dans laboratoire hydrologique d’Etat pour Eau Pota- 
ble,” Felius, LaHaye, Pays-Bas, Bulletin E., 1953, pp. 225-227. 


“Recent Ground Water Investigations the Netherlands,” W.F.J.M. Krul and 
F.A. Liefrinck, Elsevier Publishing Company, New York, 1946, pp. 

pour des problemes simultane des eaux souter- 
raine douces saluées, Santing, Commission des Eaux Souterraine, Intl, Union. 
Geod. and Geophysics, Ninth General Assembly, vol. Brussels, 1951, pp. 184-193. 

“Unsteady Flow Porous Media Means Hele-Shaw Viscous Fluid Model,” 
Todd, Trans. Amer. Geophys. Union, vol. 35, 1954, pp. 905-916. 

“Flow Porous Media Studied Hele-Shaw Channel,” Todd, Civil En- 
gineering, vol. 25, no. 1955, 85. 

“Ground Water Hydrology,” Todd, John Wiley Sons, New York, 1959, 
pp. 307-325. 

“Report Investigations Ground Water Movement and Related Problems 
Zealand Flanders Means Horizontal Slit Model (Viscous Flow Analogy),” 
Bear and Van Overstraaten Kruysse, M.P.C., M.Sc. Thesis, Israel Institute Technolo- 
gy, Israel, 1956. (Report investigations conducted the Government Institute for 
Water Supply, The Hague, Netherlands.) 

Horizontal Scale Model Based the Viscous Flow Analogy for Studying Ground 


Water Flow Aquifer Having Storage,” Santing, I.A.S.H., 
Toronto, vol. II, 1957, pp. 105-114. 
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VISCOUS ANALOGY MODELS 
DERIVATION MODEL SCALES 


The two dimensional flow ground water with phreatic surface non- 
isotropic aquifer described by: 


For point the water table, 


the capillary zone neglected and there movement across the free 
surface (See Appendix I). The same two equations can also written for the 
viscous flow the capillary interspace between two vertical parallel plates 
viscous fluid model. The only difference stems from the fact that the 


and 


substituting 


into Eq. and comparing the result with Eq. the following relations are ob- 
tained according the laws similitude: 


oa 866 


Here can any number. From the definition follows that 
therefore, the model scales are given by; 


February, 1960 
and, 
(8) 


LENGTH SCALES 


Eq. describes the distortion length scales cases nonisotropic 
When Kxp the horizontal and vertical scales must necessarily 
the same. Very often model must imitate aquifer having length much 
longer than its depth (for example, several kilometers compared hun- 
dred meters). impossible design such model without distortion 
scales. such cases, necessary assume that the hydraulic conductiv- 
ity the horizontal direction much larger than the vertical direc- 
tion, and then the distortion the scales will determined accordingly. The 
error which might thus introduced (in case disagreement between the 
assumed and the existing hydraulic conductivities) usually small. One has 
only remember that the flow essentially horizontal with relatively 
small vertical component due the elongated shape the aquifer. However, 
there are cases, such the vicinity pumping wells near the outflow 
into the sea, where this assumption might yield large errors. 


TIME SCALES 


The choice the time scale will governed the time interval between 
two successive readings discharge, levels, etc. Rewriting Eq. ina more 
detailed form, 


This indicates that appropriate time scale can obtained for given length 
scales manipulating the width the interspace, and the viscosity the 
liquid, Consideration the interrelations between length and time scales 
will enable the most practical laboratory model designed for studying 
given problem. 


SCALES 


The discharge scale obtained from Darcy’s law. Starting from the pro- 
totype flow, 


and the model flow, 


follows that 


| 
t | 
a 
Z 
| 
a 
Kxr by Zr 


Similarly, 


and, follows that 


VOLUME SCALES 


The volume scale obtained from the continuity equation 
Very often areas open water such rivers, lakes, the sea have 
imitated model. thesecases the model consists three plates (Fig. 
instead two. narrow interspace kept between two the plates while 


the third one, together with cutting the middle plate the desired shape, 
serves increase the interspace the sum the capillary interspace 


Prototype 


Open water 
Back plate 
Enlarged interspace — 


Middle plate 
Capillary interspace 


plus the thickness the middle plate. When the volume liquid leaving the 
reservoir and entering the capillary interspace importance, this en- 
larged interspace subject additional condition: 

The same volume scale must applicable both the narrow and the en- 
larged interspaces. With the notation Fig. the following equations can 
written: for the soil 


and for the open water 
Where the bar indicates values for the open water. 
Since 


Vr = Vr> 


or, 


FIG. 

a 


REPLENISHMENT SCALES 
The natural replenishment scale can derived from the equation which 


approximately describes the flow phreatic aquifer with horizontal bot- 
tom. For the case shown Fig. 


or, with constant, 


where the thickness the flow and the vertical recharge rate. Since 


Kxr 


NONHOMOGENEOUS AQUIFERS 


case nonhomogeneous aquifer, that is, aquifer which consists 
several layers with different permeabilities, animitation the model pos- 
sible varying the width the interspace. can done introducing thin 
sheets plastic any other suitable material with the geometric shape 
the layers into the interspace, thereby leaving narrower space free for the 
flow the viscous liquid. 

Since the discharge scale must remain the same for the various layers 
follows that must remain constant. Hence 


might mentioned here that sometimes the condition Kxp used 
instead Eq. 15. This follows from the assumption that K,, must remain 
constant throughout the model whereas the right thing keep con- 
stant throughout the layered model. 

condition which causes some trouble when computing scales for model 
imitating layered aquifer the condition that Np/ should kept con- 
stant throughout the model. This condition follows from Eq. when kept 
constant throughout the model. Very often this condition cannot fulfilled, 
but the error introduced generally small. most practical cases, the most 
important layers with respect hydraulic conductivity are the sand sand- 
stone layers where variations permeability porosity from one layer 
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the other are relatively small that the error introduced neglecting 
this last condition small. However, larger error will introduced the 
case more impervious layers. 

the case impervious layer, therefore, the entire inter- 
space filled. When layers small permeability are present the computed 
for these layers might come out very small (for example, 0.1 mm) 
and cause difficulties constructing and operating the model. These cases 
occur very often when relatively thin semi-pervious layers, which the flow 
mainly the vertical direction, are present. practical solution this 
problem would use strips which completely seal part interspace and 
leave wider interspace the remaining part shown Fig. 

The relationship between the partly sealed and the enlarged interspace 
determined the condition that the discharge through the whole length 


Back plate 


Front plate 


the case the narrow interspace must equal the discharge through 
the shortened length with the widened interspace Therefore, 


from which follows that 


Experience shows that thecase relatively narrow interspaces the free 
space must increased about 15% over the computed values due 
increased wall friction. The use these strips was first introduced the 


Hydrological Laboratory the Government Institute for Water Supply, The 
Hague. 


STORAGE COEFFICIENT CONSIDERA TIONS 


this point only models which not imitate storage coefficient have 
been considered. was assumed that specific storage coefficient can 


| 


The problem imitating confined aquifer with specific storage coeffi- 


confined aquifer should be: 


neglected the case the phreatic aquifer. this case the storage coeffi- 
cient the aquifer identical the effective porosity whereas the specific 
storage coefficient which results from the elastic properties soil and water 
small that can assume Sop and Som and and they 
cannot used for the determination the time scale. 


cient, model has not yet been completely solved. Nevertheless, there 
are several possibilities some which have been successfully tried. 
From and follows that the design conditions for model 


(18) 


desired consider the specific storage coefficient well the 
porosity the case anunconfined aquifer, the following condition must also 
satisfied: 


aquifer’s storage coefficient results from the elasticity the soil and 
the water, but since the pressures the model are relatively small can- 


Phreatic table 


FIG, 


not use the elasticity the plates and the viscous fluid for this purpose. The 
the confined aquifer. 

One the several methods used for imitating the storage coefficient 
confined aquifer connect the back plate vertical pipe the cross sec- 
tion which will computed according the imitated storage coefficient 
sketched Fig. The aquifer length divided into several parts each with 
length (for example, cm). The storage coefficient each these 
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columns will represented vertical pipe with cross sectional area 
Hence 


When the specific storage coefficient, varies along the height the col- 
umn, possible connect every opening the back plate separate 
vertical pipe each having different diameter. this case the entire area 
the plates divided into rectangles which the storage coefficient each 


FIG, 


separately imitated. The same method can used order increase the 
storage coefficient model which imitates phreatic aquifer that 
more suitable time scale can obtained. this case 


TWO-FLUID FLOWS 


important problem which can investigated the model that the 
simultaneous flow two liquids. The problem salt water intrusion into 
coastal aquifers example one these problems. model two vis- 
cous liquids (distinguished different colors) can beused. Certain conditions 
must satisfied for this situation. 

Assuming the hydraulic conductivity the model the same for the areas 
occupied the two different liquids, and using the fact that the-hydraulic 
conductivity the soil inversely proportional the viscosity the water, 
then 


and 
m a 
| 
4 le} m Im 
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and 


where subscript indicates the fresh water and subscript indicates the salt 
water. 


table 


Fresh Water, 


Salt 


FIG, 


The pressure any point along the interface the same both sides 
the interface. Fig. showing coastal phreatic aquifer, point the 
equation for the prototype 


' 


or ap dp = hp 


where 


Similarly, for the model 


follows that the choice the two liquids will governed the condition 
that: 


that 


h’ 
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VISCOUS ANALOGY MODELS 
CAPILLARY CONSIDERATIONS 


The problem the capillary rise the viscous liquid which flowing 
the interspace must also considered. The water table this liquid does 
not represent phreatic surface (on which the pressure the atmospheric 

pressure), but the surface which results from the capillary rise. The height 
the capillary rise depends upon the width the interspace, the type 
liquid, and the material the plates. The phreatic table the model ob- 
tained deducting the capillary rise between the plates from the liquid sur- 

face elevation. 

ground water problems where the thickness the saturated layers 
much larger than the capillary rise, this phenomenon secondary impor- 
tance only. choosing suitable interspace, liquid, and vertical scale, 
the capillary rise the model can made correspond the capillary 
rise the soil. 

There are other problems, such seepage from reservoirs and channels 
the first stages the non-steady flow, where the capillary effect (in all 
directions) primary importance. Corrections can made proper 

change the driving head the reservoir the channel. 


FIG, 

example, consider the flow from artificial replenishment reser- 

voir (two-dimensional case) shown Fig. The potential point is: 


where denotes the capillary rise and the ordinates points the water 
surface. The same expressions can written both for the model and the pro- 
totype: 
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follows that 


This means that the water surface the model reservoir must raised 
that 


When raising this level, special attention must given the fact that the 
geometric boundaries through which the seepage takes place are not altered. 
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APPENDIX 


Let F(x,y,z,t) the function which describes the phreatic surface. 
the position the same particle time when the phreatic surface 
moving, then F(x Ax, Ay, Az, =0. The change can ex- 
pressed the total derivative the form: 


The potential for the phreatic surface takes the form 


Hence, 


and directions, respectively. 


According Darcy’s laws 


hence 
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where the effective porosity and all subscripts indicate directions. 
inserting (3) into (1), the resulting final equation 
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pipe area; 
width interspace; 
width investigated area prototype; 
acceletation gravity; 
capillary rise; 
hydraulic conductivity; 
hydraulic conductivity the x-direction; 
hydraulic conductivity the z-direction; 
subscript, indicates model values; 
vertical recharge rate; 
pressure; 
subscript, indicates prototype values; 
discharge; 
subscript, indicates the model-prototype ratio; 
specific storage coefficient defined the volume water released 
from storage unit volume aquifer unit decline head; di- 
mension 
storage coefficient, defined the volume water released from 
column aquifer (the height the whole aquifer) having unit cross- 
sectional area unit decline head. dimensionless; 
time; 
volume; 
constant; 
specific weight salt and fresh water, respectively; 
ratio the difference between the specific weights salt and fresh 
water the specific weight fresh water; 
kinematic viscosity; and 
ground water potential every point. 
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BOUNDARY LAYER STIMULATION RECTANGULAR CONDUITS 


SYNOPSIS 


The effects artificial stimulation the turbulent boundary layer 


were studied 0.283-ft 0.5-ft plastic conduit. 


The length con- 


duit required for fully developed pipe flow was reduced from 18.5-ft 


14-ft 1/8-in. cubes cemented the conduit entrance curves. 


INTRODUCTION 


This paper describes investigation the effects artificial stimulation 
the turbulent boundary layer rectangular conduits and indicates the ulti- 
mate use boundary layer stimulation models hydraulic structures. 

review the problem resulting the study interest. For the past 
seven years the Waterways Experiment Station has cooperated with the United 
States Army Engineer Districts field tests large hydraulic structures, 
number which have been made Sluices and conduits dams. 
The elements typical conduit, consisting the curved entrance, the con- 
trol gate section, and the conduit proper, are shown Fig. for concrete 
gravity dam. Conduits may also contain transitions downstream from the gate 
section, vertical and horizontal bends, and exit portal constrictions. flow 

passes through each section, loses energy. For convenience, the energy 

losses occurring upstream from the uniform conduit section are grouped to- 

gether and termed the intake loss. This loss considered the total 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. February, 

Analysis Section, Hydr. Analysis Branch, Hydr. Div., U.S. Army Engr., Water- 
ways Experiment Sta., Vicksburg, Miss. 

Capt., Corps Engrs., Army, Germany. 
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1.—ELEMENT TYPICAL CONDUIT 
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available energy head minus the conduit velocity head and the friction loss 
and miscellaneous losses the uniform conduit section (Fig. 2). The in- 
take loss, defined, includes the losses attributable the entrance curves, 
the gate section and slots, and the tunnel transition. may expressed 


head (conduit proper). every case where comparison between model and 
prototype results has been possible, the prototype loss coefficient 
ably greater than that observed the model. For example, the observed Pine 
Flat Dam sluice prototype coefficient 0.16. comparable laboratory study 
indicates coefficient 0.07. Similar results were noted for Denison Dam 
flood control conduits where coefficients 0.12 the model and 0.19 the 
prototype were observed. 

Models hydraulic structures are usually designed the basis geo- 
metric similitude and Froude’s law. The latter requires that the inertial and 
gravitational forces have the same relation the model the prototype. 
Consequently, the model will have smaller Reynolds number than the proto- 
type and relatively greater friction losses. Therefore, the model made 
smooth possible obtain closer friction relation the prototype. 
believed that the turbulent boundary layer may develop more rapidly the 
very high Reynolds numbers the prototype than the low Reynolds 
numbers that occur the model. Further, the entrance, gate slots and tran- 
sition the prototype couldestablish higher turbulent energy level than could 
sustained wall roughness alone. the case the model, the boundary 
layer may develop more slowly and may not reach turbulent energy level 
comparable that the prototype. such the case, artificial stimulation 
turbulence models may necessary reproduce the intake loss coef- 
ficient the prototype. 

understand the reasons forthe apparent differences between 
prototype intake losses, the Waterways Experiment Station has under- 
taken three approaches the problem. First, theoretical studies boundary 
layer development and wall shear stress curved surfaces were initiated, 
and are still progress. Second, cooperation with the Omaha District, 
horizontal strut was installed one the 22-ft-diameter flood control tunnels 
Ft. Randall Dam. This strut contained numerous instruments for measuring 
velocity distribution and turbulence fluctuations. Tests were made with flow 
fps. addition, specially designed probe was used investigate 
velocity distribution within one foot the tunnel wall. And third, laboratory 
study the effects artificial stimulation the turbulent boundary layer 


flow rectangular conduit was made. The latter study forms the basis for 
this paper. 


PREVIOUS STUDIES ARTIFICIAL STIMULATION 


For number years, engineers engaged intesting the drag onship models 
the David Taylor Model Basin, Navy, have artificially stimulated the 
boundary layer compensate for the relative difference boundary layer 
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growth the model and prototype. Breslin and report that without 
stimulation, this difference results inappreciable error the predicted power 
requirements the prototype design speeds. This error has been attribut- 


the fact that the slower development the turbulent boundary layer 

the low model Reynolds number results smaller shearing stresses than does 

the more rapid development the high prototype Reynolds number. Prior 
was but little published literature describing attempts investigate 

the effects artificial stimulation the boundary layer itself. The work 
Breslin and Macovsky was pioneer study this field. 


THEORY TURBULENT BOUNDARY LAYER 
The theory the turbulent boundary layer has been presented numerous 


publications. Much work has been done boundary layer development flat 
plates and circular conduits, and limited open channels large 


FIG, 3.—BOUNDARY LAYER DEVELOPMENT 


width-depth ratios. extensive search published literature revealed 
record any studies conducted rectangular conduits varying ratios 
height width. 
brief summary the theory boundary layer development conduit 
pipes the flow laminar with Reynolds numbers less than 2,000; that tran- 
sition range occurs between 2,000 and 3,000; and that above Reynolds number 
3,000 the flow turbulent. 

Boundary Layer Parameters.—Turbulent flow conduit develops 
progressive thickening the boundary layer the region close the wall 
until this thickness equal tothe radius the conduit which point complete 
turbulence established. the case the rectangular conduit this thickness 


“Effects Turbulence Stimulators the Boundary Layer and Resistance Ship 
Model Detected Hot Wires,” Breslin and Macovsky, without stimu- 
lation, this difference results appreciable error the Navy Department, David Tay- 
lor Model Basin Report 724, August 1950. 
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RECTANGULAR CONDUITS 


would one-half the minimum conduit dimension. This phenomenon, for 
two-dimensional flow, illustrated Fig. which the velocity 
the undisturbed central core; represents the longitudinal component the 
velocity within the turbulent boundary layer;and denotes the thickness the 
turbulent boundary layer. 

defines the nominal thickness the turbulent boundary layer 
the thickness the zone which the major portion the velocity variation 
takes place. Velocity distribution the nominal boundary layer illustrated 
Fig. 

Two important the study the turbulent boundary layer 
are the displacement and the momentum thickness The dis- 
placement thickness usually defined the distance which the streamlines 


BOUNDARY 


FIG, 4.—VELOCITY DISTRIBUTION THE 
BOUNDARY LAYER 


have been displaced from the boundary result the turbulent boundary 
layer. The displacement can expressed as: 


Aclearer concept this parameter has been presented Lighthill5 who states 


between the surface and any streamline just outside the boundary 
layer, there must constant volume flow per unit span. This will 
the flow reduction inside the layer compensated for out- 
ward displacement such streamline through distance (which 


“Engineering Hydraulics,” Rouse, John Wiley and Sons, Inc., Y., 1950. 


Displacement Thickness,” Lighthill, Cambridge University Press, 
Journal Fluid Mechanics, Vol. Part August 1948. 


4 
—_ 
—=+=0.99 
uj 


produces flow increase since the velocity the region 
streamline displacement).” 


the foregoing, the displacement thickness and the streamline ve- 
locity. 

The momentum thickness generally defined measure the defect 
the momentum flux decrease energy within the boundary layer. The 
momentum thickness can expressed as: 


Fig. illustrates the physical and mathematical significance these two pa- 
rameters. 

Other useful parameters the study the boundary layer are the local 
momentum Reynolds number: 


inwhich the kinematic viscosity the fluid; and the local resistance coef- 
ficient defined Ross:6 


DESCRIPTION TEST APPARATUS 


Tests were conducted rectangular plastic conduit 0.283 wide, 0.5 
high and long. The axis the conduit was mounted perpendicular the 
vertical face pressure tank. The transition from the pressure tank the 
conduit was made four elliptical surfaces formed the equation: 


which equal the dimension the conduit the direction concerned. 

Flow was controlled regulating valves and measured venturi meters 
located between the pressure tank and the water supply system. The pressure 
tank head was measured mercurial U-tube manometer and referenced 
the center line the conduit. The conduit pressure gradient was determined 
side center-line piezometers the conduit proper connected manome- 
ter bank. Fig. schematic plan the general test installation shown 
Fig. 

Horizontal traverses selected stations were made means stag- 
nation tube. This tube consisted traversing rod 0.25-in. diameter 
which was fixed 0.5 in.-long hypodermic needle, 0.028-in. O.D. and 0.016 in. 
The traversing rod and hypodermic needle were inserted the conduit 
the selected station prior beginning operation. Leakage around the 
traversing rod was controlled means packing glands. Fig. shows the 
traversing rod and hypodermic needle inserted the conduit. 


Flow the Entrance Region Pipe, Ross, Paper No. 
ASME Preprint, 
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FIG, 6.—PLAN TEST INSTALLATION 
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8.—STAGNATION TUBE 
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The boundary layer was stimulated blocks cemented the curved en- 
trance shown Fig. The blocks were 0.02-ft square with 0.02-ft clear- 
ance between each block. Separate made 1/16-in. and 1/8-in. 
high blocks. each test the stimulators were fastened all four surfaces 
the entrance curves. 


TEST PROCEDURE 


The procedure followed was identical for each the three series tests 
(no stimulation; 1/16-in. stimulation, and with 1/8-in. stimulation). Tests were 


ap 
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FIG, 9.—ARTIFICIAL STIMULATION 


made with heads ftand ftof water the pressure tank for each degree 
stimulation. The resulting discharge through the conduit was determined. 
The conduit side center-line piezometric pressures were read the manome- 
ter bank the nearest 0.005 ft. Velocity distribution data the boundary 
layer were obtained traversing the stagnation tube across the conduit. The 
first two readings were taken 0.005-ft intervals from the conduit wall. Suc- 
cessive readings were made 0.01-ft intervals the centerline the con- 
duit. The tube was always oriented permit determination the x-velocity 


component. Fig. shows the piezometer locations. The velocity traverse 
stations are shown Fig. 11. 
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DISCUSSION DATA 


the interest brevity, the discussion data limited tests flows 


having Reynolds number 6.02 105. This Reynolds number was computed 
using the diameter pipe equal times the hydraulic radius the test 
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FIG, 11.—VELOCITY TRAVERSE STATIONS 


conduit. The results testson flows with other Reynolds numbers were simi- 
lar the results presented this paper. 

Velocity Data.—The local static pressure each traversing station was 
computed from the pressure data measured the side center-line piezometers. 
The and the local static pressure used compute the 
velocity distribution across the conduit. semi-log plot the velocity data 
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RECTANGULAR CONDUITS 


with respect the conduit wall was made. Two straight lines could drawn 
through each set plotted data: sloping line starting the conduit wall and 
line perpendicular the conduit wall. The sloping line defines the logarith- 
mic law the velocity distribution the turbulent boundary layer and the 
perpendicular line the velocity the undisturbed core. The boundary layer 
thickness was determined graphically the intersection these two lines. 
Velocity traverses for unstimulated flow are shown Fig. 12. The pro- 
gressive thickening the boundary layer from station station shown 
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FIG, 12.—BOUNDARY LAYER UNSTIMULATED 
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FIG, 13.—BOUNDARY LAYER STIMULATED 


the position the intersection the velocity profile lines. The absence 
the constant velocity line station indicates that the boundary layer 
fully developed. This station located about equivalent diameters down- 
stream from the entrance. 

Velocity traverses for stimulated flow are shown Fig. 13. With 1/16-in. 
block stimulation, the sloping line station extends the center line the 


conduit. almost similar condition exists station with 1/8-in. block 
stimulation. 


The boundary layer thicknesses each station was determined from the 
plots the velocity profiles. These data, plotted Fig. 14, show the pro- 
gressive thickening the turbulent boundary layers with and without stimu- 
lation. The effects successive increases the degree stimulation are 
shown the inward movement the boundary layer and the decrease con- 
duit length required for fully developed turbulence. 

Pressure Data.—The observed conduit pressure profile without boundary 
layer stimulation shown Fig. 15(pressures referenced side center line). 
This profile indicates that the gradient extends within11 
the conduit entrance. The steeper the pressure gradient upstream 
from this point occurs the regionof the conduit where the boundary layer 
rapidly developing. The effect stimulationon the pressure gradient shown 
Fig. (pressures referenced side center line). With each increase 
the degree stimulation, corresponding decrease occurs the elevation 
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FIG, 14.—STIMULATION AND BOUNDARY THICKNESS 


the pressure gradient the regionof boundary layer development together with 

slight upstream extension the friction pressure gradient line. 
Friction.—The observed pressure gradient the last ftof the conduit was 

used compute friction factor the Darcy-Weisbach formula. For flow 

with Reynolds number 6.02 105 “f” value 0.011 was obtained. This 

value plots below the smooth pipe curve the Moody-type friction diagram. 
The equation for the wall-shear-stress coefficient, developed Ross 


for circular pipes was used show qualitative effects stimulation the 
wall shear stress. Application this equation the data from the unstimu- 
lated boundary layer resultedin decrease the wall-shear-stress coefficient 
from 0.0027 station 0.0022 station where the boundary layer was 
fully developed. With maximum boundary layer stimulation the coefficients 
were found decrease from 0.0026 station 0.0023 station again 
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FIG, 15.—PRESSURE PROFILE STIMULATION 
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where the boundary layer was assumed fully developed. While the Ross equation 
was developed from data for circular pipes, its application rectangular 
conduit considered valid. 

Intake Losses.—The empirical intake loss coefficient definedearlier inthis 
paper dependent upon extension the friction gradient upstream the in- 
take structure. slight decrease the slope and the friction pres- 
sure gradient was noted with maximum stimulation. believed that this re- 
sulted increase the intake loss and decrease discharge too small 
measure with available equipment. hoped that other tests can made 
with greater stimulation and more refined equipment. 


The following conclusions can drawn from the investigation: 


Unstimulated fully developed pipe flow will occur about equivalent 
diameters from the entrance rectangular conduits with height-width ratio 
1.77 when operated the range Reynolds number 

Stimulation the boundary layer will effectively reduce the conduit 
length required for fully developed turbulent flow about equivalent diame- 
ters. 

The friction pressure gradient indicates fully developed turbulence 
about diameters from the entrance. 

Computation intake losses based pressure gradients within di- 
ameters the entrance will result erroneous indications the intake losses. 
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SCOUR BRIDGE CROSSINGS 


SYNOPSIS 


Relationships are proposed for the prediction scour piers and 


abutments for the case which sediment supplied hole. 


The relationships were obtained from combination approximate 


analysis and laboratory experiments, and depend knowledge the 


flow conditions the bridge site. 


GENERAL ASPECTS THE PROBLEM 


order design the foundations bridge over alluvial stream, 

necessary know the elevation the stream bed the vicinity the 
piers and abutments. The elevation question, however, not the elevation 
the time survey happened made, but the lowest elevation which will 
occur during the anticipated life the bridge. Unfortunately, engineers are 
equipped with slide rules rather than crystal balls; therefore this lowest 
elevation must coupled with probability occurence. this basis 
economic analysis possible with the cost construction for one stream 
bed elevation lower than another considered insurance premium for the 
decreased chance loss (including interruption traffic). 

Causes Scour Bridge Piers and Abutments lowering the stream 
bed the vicinity the piers and abutments can occur from variety 
causes. useful distinction can made separating the various causes 

into two general categories; (1) those characteristic the stream itself, and 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. February, 1960. 
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(2) those due the modification the flow the bridge crossing. For the 
past ten years scour due the modification the flow has been studied the 
Iowa Institute Hydraulic Research under the sponsorship the Iowa State 
Highway Commission and the Bureau Public Roads. result that in- 
vestigation means for predicting the scour bridge piers have been evolved. 
The relationships which were obtained will presented skeleton form 
later section; for full discussion the investigation and the design rela- 
tionships, reference should made the final reports which have been pub- 
lished bulletin form the Highway Research Board.2,3 

Even without the complicating feature bridge crossing the bed and 
banks alluvial stream cannot considered fixed. Superposed the 
very slow changes measured geological time, there are changes which may 
occur suddenly single flood, periodically during series water years, 
slowly during the life bridge. Past observations the chosen site, 
similar sites, must relied upon estimate the scour which may occur 
the future. Each river and each reach must studied understand its in- 
dividual, almost personalized, characteristics. 

The reach which includes the chosen bridge site may degrading: 
naturally, the erosional agent the geological cycle, consequence 
dam some distance upstream, the result stream straightening 
downstream. Old surveys and other records should yield evidence natural 
degradation which can used predict the future. Degradation caused 
the works man can analyzed, least approximately, the basis 
sediment-transport relationships, that is, the plans man can antici- 
pated. 

Meandering streams will tend have the greatest depth the outside 
the bend and shallows between the bends. During high water, the bends will 
scour and the crossings fill; during low water, the roles will reverse. the 
loops the meanders grow, there may natural cutoffs with resultant deg- 
radation the reach above and aggradation the reach below. Usually the 
meander grows, scouring the outside the bend, filling the inside; not 
unusual during flood, however, for chute develop across the inside 
bend and even replace the old channel. 

Braided streams are characteristically wide and shallow, with deeper 
channels which may shift erratically. Even not believed that any- 
thing can wrong, will” natural law, must admitted that there 
fair chance that during the life bridge deep channel may develop 
pier the time flood. 

Any contraction, whether the main channel the overbank flow, will 
result scour during high stages. shown the principle 
continuity applied both the discharge and the sediment load permit this 
case solved analytically since the flow conditions both the contracted 
and uncontracted reaches can considered uniform. likely that this 
case scour has given rise the old riverman’s belief that during flood 
plains river will scour its bed much the water surface rises. That such 
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scour cannot general throughout the stream has been well demonstrated 

The Lang Contraction.—By describing the conditions the uniform 
reaches above and long contraction discharge equation and 
sediment-transport equation, possible solve for the ratio the depths 
flow the two reaches—and, therefore, the depth scour. The solution 
will depend, detail, upon the equations selected describe the flow and 
transport. Fig. definition sketch the long contraction which there 
both channel and overbank contraction. The Manning formula can used 
describe the flow conditions 
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approximate form the total-sediment-load relationship recently pro- 
posed6 can used describe the sediment concentration 


which 


That portion the shear associated with the sediment particle size 
termed and the critical shear for sediment movement. The coef- 
ficient and the exponent depend the mode sediment movement and 
are functions the ratio shear-velocity fall-velocity. Only the exponent 


importance the final solution and the functional relationship can 
expressed 
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For flood flow conditions, the assumption can made without serious error 


algebraic manipulation then permits the depth ratio expressed 
function the contraction parameters 
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FIG, 1.—DEFINITION SKETCH LONG CONTRACTION 
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FIG. 3.—EQUILIBRIUM DEPTH SCOUR PIER 


: 
2 
1.5 
1.0 


Since the ratio the values should not too different from unity and the 
power the most 0.37, this factor can safely neglected. Assuming that 


the material scoured out during flood redeposited over large area 


that the depth scour the foregoing equation for overbank 
construction reduces 


and for channel contraction 


0.59 


B,\ 0.69 


These relationships are presented graphically Fig. 

The case the long contraction interest not only for its own sake, 
but for the light its solution sheds the nature scour rivers. The 
shear-velocity the contracted and uncontracted reaches particular 
stream should not very different; that for given mode sediment 
movement the depth scour for channel contraction is, first approxi- 
mation, dependent only the geometry the contraction described the 
width ratio and the approach depth. The assumptions which give rise the 
approximation are that the value the same the two reaches, and that the 
ratio large compared unity. Although reasonable and justifiable, 
these assumptions result the neglect any secondary effect velocity 
and sediment size. 

the case the overbank constriction, the ratio the total discharge 
the channel discharge and the depth the approach flow are sufficient de- 
termine the depth scour, granted the two assumptions previously mentioned. 
this case, the mode sediment movement does not affect the depth 
scour. The detailed geometry the approach conditions will, course, de- 
termine the flow distribution determining the velocity and cross-sectional 
area the flow the channel and the floodplain. However, matters not 
whether the velocities are high low, whether the flow sections are large 
small; but only what discharge ratio results therefrom. 
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bridge crossing effect long contraction foreshortened such 
extreme that has only beginning and end. The flow the crossing can- 
not considered uniform, but the solutions for the long contraction can 
modified describe the scour bridge piers and abutments with the use 
experimentally determined coefficients. 


LOCAL SCOUR PIERS AND ABUTMENTS 


Experimental Investigation.—Following the floods 1947 which the 
State experienced considerable monetary loss because thedam- 
age the failure bridges, the Iowa Institute Hydraulic Research em- 
barked investigation scour around bridge piers and abutments. For 
the preliminary phase the study which the effect the geometry the 
pier abutment was studied, flume 10.5 wide and long was con- 
structed. Two supply lines and two tailgates permitted the wide flume 
operated two flumes the addition center wall. Since neither 
sand elevator nor trap was built into the flume, the time period the runs 
was limited and only qualitative results could obtained. The velocity and 
depth flow the standard run, however, resulted general movement 
the bed, and the observations which were made established the nature the 
phenomenon that was occurring. 

Albeit qualitatively, the importance the length-width ratio the piers, 
the angle attack the stream against the piers, and the length encroach- 
ment the abutment was established this flume, well the unimpor- 
tance small details the geometry pier abutment. was also dem- 
onstrated experimentally that there was equilibrium limiting depth 
scour. the preliminary work, piers and abutments representative cur- 
rent Iowa designs were studied. The study geometry was continued? with 
more generalized pier forms and scour arrestors. The effect debris was 
was also investigated (debris can considered enlargement the pier 
somewhat less than the size the debris mass depending the permeability 
the mass). 

For the second phase the program which the effects velocity and 
depth flow and sediment size were assessed, transport flume was 
constructed which was wide and long and equipped with sand ele- 
vator and trap. Most the runs were made with representative Iowa pier 
set angle 30° the flow. The depth scour was measured 
electrical scour meter which could sense the position the sand the dif- 
ference conductivity between the water-sand mixture and plain water. 
wide range velocity, depth, and sediment size was utilized the charac- 
teristics the flume permitted. Froude numbers approaching the critical 
and sand jumping the trap limited the maximum value the velocity and the 
minimum size the sand, respectively. The minimum value velocity and 
the maximum size the sand were limited the requirement that there 
general movement the sediment bed load. Nevertheless, sixty-fold 
change the rate sediment transport was obtained—probably the most sig- 
nificant indicator the range the tests. The depth scour the front 
the pier shown Fig. systematic scatter could detected the 
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that the depth scour was uniquely determined 
the geometry. 

This transport flume was also utilized for model study bridge pier 
which field measurements were made. The measurements the field were 
obtained with adaptation the laboratory scour meter. Comparison the 
model and prototype data indicated that the depth scour could treated 
simply another length and that the equilibrium depth scour obtained the 
field. 

Measurements system multiple cylinders were also made. The 
most significant findings these tests were that the depth scour did not 
depend the degree contraction (or proximity adjacent cylinder) 
until the scour holes overlapped, and that the minimum depth scour when 
interference did occur was given the solution for the long contraction. 

For the third phase the study, the first flume was modified being 
lengthened ft, and equipped with sand elevator and traps. The sand was 
supplied five one side the flume, simulating the river 
Clear water was supplied the other side the flume which simulated the 
floodplain, overbank area. The experimental program for the investigation 
the scour the abutment included variation the ratio discharge 
overbank and channel, the depth flow, the position and angle the pier, and 
the presence vegetal screening along the bank line. few experiments 
were also conducted the effect spur dikes off the end the abutment 
parallel the bank, and screening the floodplain the approach the 
abutment. 

Design Relationships For Abutments.—The observation previously men- 
tioned, that the local depth scour does not depend the degree contrac- 
tion until scour holes around neighboring obstructions overlap, suggests 
extension the solution for the long contraction the case local scour. 
For normal sand, the width scour hole right angles the flow 
about 2.75 times the depth scour. Fig. definition sketch long 
contraction, which will assumed approximates the case the abutment. 
The scour the long contraction assumed fraction 1/r the scour 
the abutment Rewritten for this special case Eq. becomes 


Since the depth scour unknown, unknown. Therefore, the discharge 

width should approximately equal 2.75 and practice trial-and- 


error procedure will necessary). Algebraic manipulation will now result 


7/6 


which the average depth flow the width Eq. 10, with value 
4.1, Fig. together with the experimental data for runs 
with the approach fill and abutment normal the direction flow. Various 
set-back distances the abutment and clear distances between the abutment 
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FIG, 5.—SCOUR OVERBANK CONSTRICTION 
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and the vegetal screen are represented among the data. The solid points are 
from the last, and best, series runs which the abutments was not set back 
from the bank line. The high points are from early runs which the velocity 
the approach the overbank was great unrealistic. The floor 
the overbank area was later lowered smaller velocity ap- 
proach. 

For large values the parameter apparent that the ap- 
proximate solution describes the scour remarkably well. For small values 
the parameter, the discrepancy due cross flow from the channel 
the overbank the approach the constriction. This cross flow could 
observed means dye streaks. Since the water-surface slope the 
channel was greater than the overbank area, the water surface the head 
end the flume was higher the channel than the overbank. Correcting 
the overbank discharge approximate calculation the cross flow would 
shift the experimental points over the analytic curve. This has not been 
done, emphasize the importance fully determining the overbank flow. 
practice, the rate flow the floodplain will determined considering 
the slope the water surface, the depth flow and the area the flow sec- 
tion, and estimating value. the overbank flow small, cross flow 
the area immediately upstream the constriction can give rise condition 
similar that the laboratory flume. this cross flow possible, but can- 
not evaluated, the upper curve Fig. should used. 

similar approximate solution can obtained for the case en- 


croaching abutment. Using Eq. for the condition bed-load movement 


1.70 


where the effective length the abutment such that (denoting the flow 
obstructed the abutment 


with equal 11.5 plotted Fig. together with some experi- 
mental points. The square points were during investigation the 
relief bridge problem and are for the case vertical, blunt-ended, normal 
wall. The round points are data from the investigation multiple cylinders 
which there was interference effect. 

Among the variations geometry tested was the angle incidence be- 
tween the direction the approach fill and the direction the flow the 
channel. The results are presentedin Fig. multiplying factor, co- 
efficient, applied the depth scour calculated for the normal 
crossing 90°). Although only determined for the overbank case, these 
coefficients should also applicable the encroaching abutment. check 
the points for the cylinders plotted Fig. would indicate that the effec- 
tive angle the curved wall the cylinder about 45°. 

The effect setting the abutment back from the normal bankof the stream 
difficult assess. the laboratory experiments measurable effect 


could noted. The long contraction solution can further modified con- 
sider that part the overbank flow remains the floodplain, and that the 
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6.—SCOUR ENCROACHING ABUTMENT 
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FIG. 8.—EFFECT THE RATIO SHEAR-VELOCITY FALL-VELOCITY 
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scour hole moves back with the abutment. This approximate solution also in- 
dicates that the effect set back small. Both the experiments and the 
analysis, however, not model real set back certain respects. Because 
the bank, the depth flow the vicinity the abutment will less. 
the set back very large, the presence the river nearby immaterial. 
The condition then special case relief bridge and the scour due 
clear water. the bank material contains considerable clay well 
covered with good, well-rooted, tight turf the material near the abutment 
might very resistant scour. 

For the case the encroaching abutment, the mode sediment movement 
will affect the depth scour. approximate evaluation this effect can 
obtained comparing Eq. 8a, andc. For the same values the 
depth scour for the different modes movement were calculated. Defin- 
ing ratio the depth scour under suspended-load conditions 


depth scour under bed-load conditions, Fig. was prepared. Note that the 
fall velocity should that the bed material being scoured out. The so- 
called wash load should ignored. The effect change the mode 
movement should applicable also the local scour pier. Experiments 
made with circular cylinder and bed material 0.04 diameter con- 
firmed the order magnitude the effect. The depth scour—very difficult 
measure and extremely sensitive unsteadiness the about 
50% greater than for bed-load conditions. The ratio shear-velocity fall- 
velocity this experiment was about 20, and the concentration between and 
10% weight, that is, 50,000 and 100,000 parts per million. The ratio be- 
tween the rate suspended-load movement and the bed-load movement was 
about 500 

Design Relationships For the available data piers were ad- 
justed scour around rectangular pier aligned with the flow, and the design 
curve shown Fig. was drawn with conservatism and with due regard for 
the reliability the various data. Also shown Fig. Eq. 
11.5 (note that The most important aspect the geometry the pier 
was the angle attack between the pier and the flow, coupled with the length- 
width ratio the pier. family curves shown Fig. multiply- 
ing factor applied the depth scour obtained from the basic 
curve. the pier set angle the flow, recommended that al- 
lowance made for shape. certain that the pier aligned with the 
flow and will remain during the life the bridge, shape factor can 
applied shown Table 

All the experimental work which the design curve was based involved 
bed-load movement. the mode movement different, the ratio shear- 
velocity fall-velocity factor Fig. should used. 


Besides the local scour the pier, possible that the scour centering 
the abutment may extend its area influence out piers near the bank. 
Thus, fraction the depth scour the abutment may have added 
the local depth scour. Typical scour holes have been shown whereby this 
effect can estimated.3 

the bulletins the Iowa Highway Research Board the application de- 
sign fully discussed, including schematic design examples illustrating the 
use all the design relationships. the basis these relationships the 
prediction local scour piers and abutments simple and straightfor- 
ward. order use the relationships, however, the flow conditions the 
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Shape coefficients for nose forms. 
(To used only for piers aligned with flow.) 


Nose form Length-width Ratio 


Rectangular 1.00 


Semicircular 0.90 


Elliptic 


Lenticular 
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bridge crossing must known. any civil-hydraulic engineer knows, the 
forecasting the river flow the future not easy. Few gaging sites have 
been operated for 100 years, and 100 years record are not sufficient de- 
fine precisely the magnitude 100-year flood. Unfortunately, even rarer 
floods during the anticipated life, say years, bridge, have probability 
occurence which large enough reckoned with. This problem, 
course, simply the old basic problem hydrology. 

addition magnitude, other characteristics the flow expected 
must foreseen; the stage, the direction flow (which may vary with stage), 
the division flow between channel and overbank (which will vary with stage 
and the stage the floodplain). Although these details can never foretold 
with complete confidence, they are essential the prediction the scour 
which will occur. The evaluation the flow characteristics must judicious, 
and calculated risk must accepted; otherwise, the worst conditions must 
assumed and overdesign will usually result. sufficient control can 
exercised over the valley upstream from the bridge crossing, the risk can 
reduced the chance the flood occurrences. Each bridge crossing will 
have treated individually, and the final analysis economic considera- 
tions will dictate the solution. 

Backwater Bridge Crossings.—When scour occurs bridge crossing, 
the flow pattern which then obtains does not resemble that two- 
dimensional constriction. The flow approaching the obstruction (pier abut- 
ment and embankment) ducks beneath the surface and passes through the con- 
striction spiral roller. The remainder the flow passes over this roller 
and through the constriction barely noticing the piers and abutments, that 
slot; open-channel constriction with inerodible boundaries. 
result, the velocity increases only must transport the sediment load 
downstream from the bridge; the energy losses and due the bridge 
crossing are minimal. Since the flow disturbed, there are, course, some 
losses. Fig. shows the measured rise water surface the abutment 
flume. This measured rise partly due the velocity-head change between 
the approach flow the channel and the flow the contracted section down- 
stream from the bridge. The average depth the downstream section can 
obtained from Eq. and the curve for the velocity-head differential has been 
evaluated this basis. The backwater the bridge the difference be- 
tween the velocity-head differential and the measured rise the head the 
flume. The backwater comparable inerodible flume according Liu, 
Bradley, and also plotted for comparison. readily apparent 
that the increase the flow section due the scour reduces the backwater 
fraction the amount that otherwise would occur. 

likely that major cause backwater bridge crossing al- 
luvial stream might the resistance flow the contracted stream the 
reach below the bridge. Depending bank height and cover the floodplain 
the cross flow from the stream back the floodplain might take place very 
slowly. the resistance flow greater the contracted reach than the 
normal stream, backwater will result. Moreover, any backwater which does 
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occur alluvial river tends cause aggradation upstream, and con- 
comitant lengthening the backwater curve. 


SUGGESTIONS FOR FURTHER STUDIES 


Field Measurements.—In order apply with confidence the relationships 
proposed for predicting the scour bridge piers and abutments, field meas- 
urements are needed verify the conformity model and prototype. Some 
these measurements should made sites extremely simple geometry 
like that the Skunk River pier near Ames, Iowa. Such sites can easily 
modeled the laboratory, and questions should arise extraneous ef- 
fects. Other measurements should made sites complex geometry 
which perhaps cannot modeled the laboratory. The scour can pre- 
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inerodible channel 


Measured water- 
Backwater 


11.—BACKWATER ALLUVIAL CHANNEL 


dicted with the aid the relationships presented, and any discrepancies may 
indicate factors that should receive further study. 

The experience the past can used substitute for field measure- 
ments. hindcasting, the scour which occurred around piers and abutments 
selected crossing during major floods can evaluated. bridges that 
failed and that did not fail during the same flood are checked, correlation 
should possible among hindcasted depth scour, permissible depth 
scour (based foundation design), and failure non-failure. 

Laboratory Investigations.—The one general aspect scour bridge 
crossings which has not been discussed the scour relief bridges. Relief 
bridges are openings placed the highway embankment crosses the 
floodplain the valley. For design purposes the flow through the relief bridge 
must considered consist clear water. there sediment load 
will tend much finer size than the material which will scoured out. 
The depth scour due clear water most assuredly function the 
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velocity flow and the sediment size well the geometry the opening 
and its environs. The limiting condition for this type scour boundary 
shear that equal the critical for the material which the boundary 
composed. Since there supply sediment coming into the scour hole, 
the limit there can sediment going out the scour hole. This limit 
finite, but reached asymptotically infinite time. very difficult 
the laboratory scale all factors controlling clear water scour and combi- 
nation experiment and analysis indicated. 

For abutment which set back from the normal bank less than the 
depth scour, the relationships proposed herein should applicable. 
abutment set back over three times the depth scour should considered 
relief bridge. Following successful investigation the relief bridge prob- 
lem, would seem order study the effect set backs between these 
limits. 

One could, course, studying the effect various geometries 
piers, abutments and overall crossing. However, the results this investiga- 
tion indicate that the effects these details are minor, and less importance 
than the error that can expected the evaluation the flow conditions 
the site. Further studies geometry should probably deferred until field 
measurements indicate they are sufficient importance. 
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SEDIMENT TRANSPORT AND DELTA FORMATION 


SYNOPSIS 
This paper deals with sedimentation computations that were carried 
out connection with dike design for the Saskatchewan Delta Reclama- 


tion Project. 


The first part the paper gives description the re- 
gime the Saskatchewan River general and the Saskatchewan Delta 


particular. 


The second part presents some observations delta 
formation. The third part discusses the sedimentation problems that 


were encountered the design dikes. 


The fourth and final part pre- 


sents the method that was adopted determine the future rate sedi- 


mentation the floodplains adjacent the projected dikes. 


REGIME SASKATCHEWAN RIVER 


The Saskatchewan River system, Fig.1, originates the Rocky Mountains, 
where some the peaks exceed elevation 10,000 ft. above sea level and 
where the permanent creeks begin flow elevation 5,000 ft. 7,000 ft. 
Most the runoff the mountains comes from the annual snowmelt June, 
which may last for several weeks. During normal snowmelt, there 
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sediment transport speak of. The water inthe streams clear andthe 
bles and boulders the stream bed not move. Most the sediment trans- 
port takes place during short, concentrated rainstorms. such time, the 
creeks become torrents and they bring down the debris from the mountain 
slopes tremendous quantities. The majority the mountain streams flow 
V-shaped valleys with steep gradient. They cascade over rapids and water- 
falls and are apparently state active erosion. Some the streams, 
however, flow valleys with mild gradient. These streams have braided 
character. The valley floor filled with gravel and small boulders and gives 
the impression being horizontal plane between the mountain sides, over 
which the numerous stream channels wander back and forth will. Ofcourse, 
this valley plane not quite horizontal, but has longitudinal slope which 
sufficient carry discharge and sediment load. 

Numerous small headwater streams leave the foothills the mountains 
elevation about 4,000 ft. Their bed consists boulders and gravel; their 
slope about ft. per during normal snowmelt conditions, the wa- 
ter has visibility from ft. ft. 

this place, one would normally expect the beginning the alluvial fans 
that are often found the foot mountain ranges. However, the Saskatche- 
wan Basin, appears that the slope the land east the mountains steep 
enough carry the streams with their entire sediment load that deposi- 
tion takes place. fact, the slope steep that active erosion takes place. 
Near the foothills, the stream valleys are but easterly direction, 
they increase depth and width and unite into valleys still greater dimen- 
tions. The visibility the water rapidly decreases from practically clear 
visibility about four ins. over distance about miles, indicating that 
sediment derived from local river-valley erosion. 

Between Macleod and distance about miles river, the 
following observations were made during normal snowmelt flood stage 
June. The river valley had width about mile, places narrowing down 
tion from about ft. 125 ft., and occasionally 175 ft. The river mean- 
dered from one valley bank another, and where touched, huge active land- 
slides could seen. The valley banks consisted shale overlain glacial 
There were signs ancient alluvial fan deposits. The river hada 
varying width several hundred feet and varying maximum depth from ft. 
ft. The bottom consisted coarse gravel and boulders, well-sorted and 
well-rounded, indicating that occasional transport this bed material takes 
place. While the boat drifted the stream, the clear tinkling sound coarse 
sand and fine gravel moving over the bed could heard, the bottom the boat 
acting soundboard. The slope the river was about per mile. The 
suspended-sediment concentration was found range from 400 ppm ppm. 
The low floodplain banks the river channel were composed layers sand 
coarse gravel, each layer well-sorted. one place, old oxbow was ob- 
served higher elevation than the present river channel, indicating that the 
river state degradation. Except for few minor rapids, the river 
had rather uniform maximum surface velocity about fps. This uniform 
may due the fact that the river valley excavated soft shale, 
the coarse bed material being transported from upstream. Thus, there 
cause for major rapids waterfalls. The gradient the river probably 
governed the one hand the amount bed-material load that comes from 
the mountains, and the other hand the rate erosion that takes place 
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the valley. would seem that the eroding shale banks form important 
source the suspended-sediment load that found the river. 

Between Lethbridge and Medicine Hat, the river valley similar the 
foregoing description. The river channel and river valley gradually become 
wider and deeper, but the bed material remains composed coarse gravel 
and small boulders, and the actively eroding valley banks are still composed 
shale with overburden glacial drift. 

Between Medicine Hat and Empress, interesting change takes place. The 
river banks increase height until they are more than 400 ft. above the river 
level. There are floodplains, and immediately adjacent the river channel, 
the valley banks, composed soft sandstone shale and glacial drift, 
rise weird formations that resemble badland topography. The slope 
the river decreases from about ft. per mile ft. per mile, and the maxi- 
mum surface velocity from fps fps. The maximum depth ranges from 
ft. ft., and the width the river from 700 ft. 800 ft. The interesting 
feature this stretch the rather abrupt change inthe composition the bed 
material. Near Medicine Hat, the river bed composed pebbles and small 
boulders. Halfway between Medicine Hat and Empress, the river bed com- 
posed small boulders, with intermittent stretches sand and fine gravel. 
Near Empress, the river bed composed sand, with occasional patch 
gravel boulder sticking out the sand. These last observations were 
made dragging sampler over the river bottom from drifting boat, each 
time over length several hundred feet. The type noise that the sampler 
made, plus the contents the sampler, provided reliable indication the 
composition the river bottom. 

Between Empress and The Forks, the river valley from one several 
miles wide and from one several hundred ft. deep. There are still flood- 
plains speak of. The rolling valley sides, covered with grass and brush- 
wood, come down the river edge. The river channel has slope ft. 
per mile and from one several thousand ft. wide. Within the rigid bound- 
aries the valley banks, the river tends braid. Even with small boat and 
during the snowmelt flood, one frequently gets stuck shallow water. The bed 
material consists rather uniform sand with mean diameter 0.30 mm. 
There are signs that erosion very active this stretch the river. 

the foregoing discussion, the course the Old Man River and the South 
Saskatchewan River was followed. The North Saskatchewan River has very 
similar characteristics. The Bow River carries less sediment, probably due 
its milder gradient and lesser valley bank erosion. The Red Deer River has 
its drainage area mostly the rolling plains east the Rocky Mountains. 
Normally, the runoff from these plains low, but during heavy rainstorms, the 
river may carry considerable discharge with relatively high sediment con- 
centration. This sediment partly produced the plains and partly the 
river valley, which has its middle part very pronounced badland topography. 

Between The Forks and the delta, the gradient the river increases 
several feet per mile. Occasional rapids are found, and the bed material con- 
sists gravel and boulders. The valley banks are from one few hundred 
feet high and are composed glacial drift. places, where the river chan- 
nel nudges into the valley banks, huge active bank slides over the total valley 
height, and several thousand feet long, can seen. the end this stretch, 
the river channel descends about ft. over Tobin and Squaw Rapids and 
ters the Saskatchewan Delta, shown Fig. 

About 10,000 yr. ago, receding glaciers allowed the Saskatchewan River 
take its present course. Near its lower end, where the river entered glacial 
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Lake Agassiz, most the sediment was deposited and the formation the del- 
area, shown the present time, was commenced. Deposits 
coarse material were confined channel locations, while the fine silts and 
clays settled down farther away from the river. After the channels were ex- 
tended over some distance into the lake, while possibly being split into 
several branches, situation developed whereby the river would somewhere 
break through its own formed banks and choose entirely new course. Inthis 
fashion, the Saskatchewan River changed from one group another 
for thousands years. Some the oldest channels are completely covered 
subsequent sedimentation. Some the more recently abandoned channels are 
still plainly visible aerial photographs and show the landscape 
densely wooded ridges, enclosing low-lying marsh areas. The whole delta 
area, disregarding local ridges and depressions, gently sloping plain 
about miles wide and 120 miles long, with gradient about ft. per mile 
near the upper end and 1/2 ft. per mile near the lower end. The single river 
channel that enters the delta cone has slope about ft. per mile and may 
carry during normal flood stages peak flow 100,000 cfs, containing 0.3% 
weight sediment. The several river channels that leave the delta and flow 
into Cedar Lake, one the remainders the ancient glacial lake, have slope 
about 1/4 ft. per mile and may carry during normal flood stages combined 
peak flow 50,000 cfs, containing 0.1% sediment. The difference sediment 
load between the upper and lower end used raising the delta surface. The 
remainder the sediment load used river-mouth extension into Cedar 
Lake. 

The average annual sediment transport into the Saskatchewan Delta esti- 
mated 12,000 acre-ft. About one third this volume sand with average 
diameter 0.20 mm., and the remainder wash load (silt and clay). Prac- 
tically all this sediment transported suspended state and deposited 
somewhere the delta. During low and medium flood stages, the flow wa- 
ter and sediment remains within the natural river banks and deposition takes 
place the river channel and the mouth the river channels. During high 
and extreme flood stages, the river banks are overtopped and deposition takes 
place throughout the entire delta area. These three modes sediment depo- 
sition: inthe river channels, the floodplains, and the river mouth, will 
discussed more detail the following paragraphs. 

The river channels the apex the delta cone have slope about ft. 
per mile, which accordance with the observed regime slope the lower 
North and South Saskatchewan Rivers, noted earlier. The river channels 
near the foot the delta cone Cedar Lake have slope about 1/4 ft. per 
mile. This reduction river slope corresponds with reduction sediment- 
transporting capacity. estimated from observations and supporting com- 
putations that the capacity transport the washload hardly affected, but 
that the capacity transport the sand reduced about 80%. Asa result, 
sand deposited the river channels; the coarsest fractions upstream, and 
the finer fractions downstream. This reflected the composition the 
bed material, which shows mean diameter 0.30 mm., 0.24 mm., and 0.19 
mm., respectively, near the apex, the middle, and the foot the delta. 

During high flood stages, the river banks are overtopped and flow wa- 
ter and sediment directed towards the marshes. Immediately beyond the 
edge the river bank, the coarsest sediment begins precipitate. Soil sam- 
ple analyses show that the natural levees are composed the coarsest frac- 
tion the suspended sediment that carried the river. Away from the 
river, the deposited sediment becomes finer and finer, until only the wash load 
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FIG. 2.—SASKATCHEWAN 
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remains which the quiet marshes adjacent the river 
channels. After delay, which may range from several weeks, 
the overflow water, now deprived its sediment load, returns some distance 
downstream via the outlet creeks the marshes into the river channels and 
helps reduce the average suspended-sediment concentration the river 
flow. 

Cumberland Lake and Cedar Lake, river channels emerge into 
body water. Due the inertia the moving water, the river flows forward 
homogenous body water into the lake and consequently, retains its 
silt load for while. After certain distance, however, friction overcomes 
inertia, the water slows down, and sediment precipitates. This deposition 
sediment takes place where velocities decelerate, which happens mostly 
the flanks and front the retarding river current. Thus, two longitudi- 
nal bars and number shoals are formed. From aerial photographs, taken 
yr. apart, appears that the bars and shoals keep growing elevation and 
that the openings between develop into river channels. Some the less im- 
portant channels silt until only one few main channels remain. All 
the sand and most the silt carried the river deposited near the river 
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DELTA 


mouth. Some silt and most the clay deposited further away the lake. 
The bottom Cumberland Lake and the western part Cedar Lake, for in- 
stance, are almost entirely covered with soft mud. Since Cumberland Lake 
only ft. deep flood stages, occasional winds keep stirring the bottom 
deposits and result, the water that flows out Cumberland Lake always 
discolored. Cedar Lake, however, larger and deeper, causing all sediment 


permanently settle down, that the outlet Cedar Lake has clear water 
all times. 


DELTA FORMATION 


The essential process the formation riverdelta may summarized 
follows: river brings water and sediment into lake sea. Through de- 
position the sediment, the river mouth extended. Since the river needs 
slope carry its water and sediment, the lengthening the river channel will 
result rise river stages any given point the lower reach the 
river. There limit how much river can rise above the adjacent land, 
and after this limit reached, the river breaks through its own banks, chooses 
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new course, and beginsto repeat the former procedure. The manner which 
this process takes place will more detail the following para- 
graphs. 

When river enters lake the flow water not dispersed all direc- 
tions, but due its momentum, maintains for some distance its original di- 
rection and moves well-defined stream into the lake. appreciable 
amount sediment, carried the river, will deposited the sides this 
stream, where moving and still water meet, and will form two banks that pro- 
vide further guidance the river current. When after recession the flood 
the lake and river levels drop, these banks may become exposed. This will 
enable them consolidate and possibly grow some vegetation. When the 
next flood comes, the same cause for further deposition the banks still 
present and may even increased due the vegetation, while the causes for 
erosion are lessened through the consolidation presence the vegeta- 
tive cover. When this goes from flood flood and avulsion through the 
natural river bank takes place, the river mouth will extend like finger the 
lake. Sometimes, this form river-mouth extension very pronounced, and 
other times, hardly noticeable. There are several factors which affect 
this phenomenon. First all, the type sediment carried the river im- 
portant. large proportion fine sediment will produce relatively high 
sediment concentration the overflow water and rapid sedimentation 
the banks. Moreover, fine sediment will result cohesive banks that may 
delay the occurrence avulsion and thus promote the extension single 
mouth into the lake sea. Secondly, rapid progress the river mouth will 
excessive concave upwardcurvature the streamprofile. Inother 
words, the slope nearthe mouthbecomes relatively low. mouth may 
protrude over large distance lake seabefore the height the river 
banks becomes critical for avulsion. There are several possible causes for 
rapid progress the river mouth, like shallow lake large concentra- 
tion sediment. Since the slope graded river determined the bed- 
material load and not the wash load, the circumstances would ideal for 
finger-like extension when river witha low concentration bed-material load 
and high concentration wash load would enter shallow lake. The mouth 
would progress rapidly, the banks would cohesive, and the slope would 
small. 

After avulsion takes place, the river will using boththe oldandthe new 
channel for certain length time. However, there tendency for the old 
channel silt and abandoned for the following reasons. First all, 
avulsion often takes place the outside river bend where the water 
piled during flood stages. This meansthat the new channel assumes more 
less straight course, whereas the old channel turns off the left the 
right. Due the tendency the bed load move towards the inside bend, 
the old channel will receive proportionally more sediment than water. 
der carry this sediment, the river needs steeper slope, which can obtain 
only silting its upper end. Even when the distribution water and sedi- 
ment over the old and the new channel proportion, then the new slope 
the old channel still needs steeper because two separated channels are 
less efficient transporting sediment than one combined channel. 

addition the above causes for the actual silting the old channel, 
there the tendency the new channel produce local steep gradient 
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the site the avulsion. This first all because the local difference 
elevation between the river channel and adjacent land, and secondly, because 
the distance from the avulsion the point where the sea lake level 
reached usually less the new channel, especially deposits the old 
river project far into the sea lake. This steeper gradient the new chan- 
nel near the avulsion will cause draw-down water profile the river up- 
stream. This will result local scour and further decrease river stages 
near the entrance the old channel, thus depriving this channel more and 
more its original flow. 

The deposition sediment the old channel will take place such way 
that the remaining cross section equilibrium, discussed previously. 
That means that deposition will not only take place the bottom the chan- 
nel but also the sides. The channel will become shallower and narrower. 
Finally, the entrance the channel may become narrow that plugs 
with debris may become high that becomes grown over with vege- 
tation during low-water periods. After this has occurred, the channel will 
abandoned and filled with marsh vegetation. 

The rapidity this silting-up process depends upon the proportion ex- 
cess bed-material load that enters the channel. this low, for instance, due 
the entrance being located nearthe convex side the main current, the river 
will silt slowly. the meantime, the new channel also silts due nor- 
mal aggradation. Hence, the initial degradation upstream the avulsion 
checked and reversed aggradation. may even that the aggradation 
the new channel catches withthe aggradation the oldchanneland that both 
rivers become equal partners. The aspect that these two rivers are less ef- 
ficient their sediment transport than the former one river channel causes 
their slopes steeper. Consequently, the river channel upstream the 
avulsion will aggrade above the level that had the time the avulsion. 
After this has continued for awhile, another avulsion may take place, elimi- 
nating the two old channels and replacing them one new channel. 

The formation and aggradation the new channel, which initially flows 
through marshes and shallow lakes adjacent the old channel, somewhat 
similar the formation the old channel, discussed the above. How- 
ever, one point difference interest. the old channel, was assumed 
that the river bed consisted normal river bed-material. The new channel, 
flowing through low areas where formerly fine sediment was deposited, may 
find bottom cohesive material that resists scour. Hence, the formation 
the cross section may different from what would the river flowed 
over its own bed-material deposits. could expected that the beginning 
the channel will relatively wide. After the river digs into the clay deposits 
part way, could expected that normal meandering casual shifting will 
less pronounced. could also expected that the river slope can 
steeper than would accordance with the transported bed-material load, 
due its resistance against scour. 

From the foregoing discussion, follows that avulsion may lead the 
presence two more less stable and more less active river branches. 
this repeated number times, the familiar picture river arms ina 
delta, comparable the veins leaf, obtained. downstream direction, 


the river channel will become increasingly smaller due repeated bifurcat- 
ing. 
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There another reason why river channels delta decrease size 
downstream direction. During flood stages, when the river carries most its 
sediment and when the formation the channels takes place, water consistent- 
overflows the river banks towards the adjacent marshes and lakes. Asa 
result, the flow that remains the channel will become increasingly smaller 
downstream direction. Since channel formed accordance with the 
volume its flow, follows that the channel capacity will also decrease 
downstream direction. 

The consequence this decrease channel capacity that only one flood 
flow with certain magnitude will cause uniform overflow over the delta. 
Floods greater magnitude, which occur less frequently but which have 
higher concentration sediment and are longer duration, will flood more 
severely the upper part the delta area; while floods relatively small 
magnitude will overflow more extensively the lower part the delta area. 
this way, the sedimentation kept balance. 


SEDIMENTATION PROBLEMS 


was found from soil surveys that practically all the alluvial deposits 
the Saskatchewan Delta, covering area about 2,000,000 acres, are suit- 
able for agricultural development provided the water table can controlled. 
area 150,000 acres, near The Pas, has recently been reclaimed, while 
plans are being made for the future reclamation additional 1,000,000 
acres delta land. 

Reclamation part the delta surface will result smaller area being 
available for deposition sediment and hence, faster rise the surface 
level. This turn will increase the flood elevations and thus endanger the 
projected dikes. The problem arises determine the future sedimentation 
the floodplain quantitive terms and find possible remedial measures. 

Generally speaking, there are three principal methods coping with the 
present sedimentation problem. The first method reduction the sedi- 
ment load through soil conservation. The second method would aim increas- 
ing the sediment-transporting capacity the river channel river-training 
The third method includes all plans for controlled storage sediment 
reservoirs, the delta area. the following paragraphs, the applica- 
bility these methods will discussed. 

cases where the sediment, carried the river, finds origin surface 
erosion the drainage basin, methods could effective. 
the one hand, they improve the usefulness the land, and the other hand, 
they may provide solution the sediment problem. However, the sediment 
carried the Saskatchewan River finds its origin not only onthe prairie lands, 
but also stream-valley erosion the upper part the watershed. re- 
duce this type sediment supply would call for extensive channel-protection 
works. Without entering into detail, can stated that the benefits the 
delta area would too small justify the cost such works. 

The method increasing the sediment-transporting capacity river- 
training works applied extensively rivers Europe and occasionally 
rivers North America. most cases, these works also serve the purpose 
maintaining navigation channel and protecting valuable shoreline property. 
the Saskatchewan Delta, such works, which may cost anywhere from $50,000 
$500,000 per mile river channel, would too costly justified the 
benefits, least for the time being. 
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Hence, the sediment cannot retained the place its origin, 
ried through the mouth the river, then must deposited somewhere 
along the course the river, and the problem find the location that will 
the least damage. Storage upstream sediment-control reservoir would 
provide adequate solution. For the first few hundred years, the sedimenta- 
tion problem the Saskatchewan Delta will not severe justify the 
construction exclusive sediment-control reservoir. However, upstream 
reservoirs are going built for other purposes, the feature sediment 
control will added automatically. could even expected that the conse- 
quent reduction sediment load will cause degradation the river channels 
the delta, and thenceforth, will lower flood stages and improve drainage. 

upstream reservoirs are going built during the first few hundred 
years, the Saskatchewan River will deposit its sediment load the delta. 
the present time, the larger part the sediment load deposited the Cum- 
berland Lake Delta. was from quantitative sedimentation study that this 
situation can maintained for another 150 yr. From then on, increasing 
quantities sediment will pass the Cumberland Lake Delta and will de- 
posited stream channels and floodplains downstream. For some time 
come, the only way cope with the consequent rise flood levels will 
raise the adjacent dikes accordingly. 

There undoubtedly limit the possibility raising dikes keep pace 
with the natural sedimentation outside the dikes. The ultimate slope the 
river whereby can carry the complete sediment load throughout the delta 
area about ft. per mile. line under this slope, beginning the future 
mouth the river near the outlet Cedar Lake, will reach point near 
Cumberland House more than 100 ft. above its present elevation. appears, 
therefore, that the area which presently going reclaimed will have 
sacrificed sometime the far future. However, would seem entirely 
feasible compensate for this loss reclaiming the areas that have been 
silted the meantime. Thus, situation would develop whereby every 
thousand years so, another section the delta area made available for 
reclamation. When spread over periods thousand years, the cost mov- 
ing the inhabitants from the area, which going back the river system into 
the area that going reclaimed, appears very small and should 
objection present reclamation. 

After having decided that the most appropriate remedial measure cope 
with future sedimentation the floodplain would raise the dikes need- 


ed, the remaining problem was todetermine the rate sedimentation adjacent 
the dikes. 


SEDIMENTATION COMPUTATIONS 


serious attempt was made compute detail how the future deposition 
sediment the Saskatchewan Delta would take place. Precise sediment 
sampling was carried out dozen stations. Analytical computations were 
applied determine the existing mode transport and rate 
deposition all sediment fractions for new equilibrium conditions. However, 
the end, many assumptions had made that the reliability the an- 
swer became highly problematical, particularly because was found that slight 
changes the assumptions, well within reason, could 
changes the end result. was, therefore, decided determine the rate 
future sedimentation, not the basis precise equilibrium calculation, but 
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OISTANCE IN MILES 


3.—CROSS SECTION A-A 


the broad basis estimated average annual sediment load coming into 
the delta and estimated distribution deposition over the delta surface. 
Before these two items are discussed some detail the following para- 
graphs, may appropriate point out this place that detailed analysis 
and precise calculations were definitely helpful exploring the nature the 
sediment problem, and providing insight into the proper solution the prob- 
lem. However, the present sediment problems appeared too complicated in~ 
sist precision the very last answer. 

The average annual sediment load was estimated the basis sediment 
sampling carried out during yrs. relationship was established between 
river flow and total sediment transport. From this relationship and the flow 
duration curve, sediment-transport duration curve was prepared, which 
turn yielded the average annual sediment-transport figure 8,000 acre-ft. 
per yr. Since the yr. observation were preceded about yr. above- 
normal rainfall conditions over the drainage basin, resulting better-than- 
normal vegetative cover over the sediment-producing prairies, was decided 
increase this figure 50%, resulting total average annual sediment 
transport 12,000 acre-ft. per yr. 

The average distribution sediment deposition over the delta area was 
estimated from study profiles and cross sections delta cone. One 
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such cross section shown Fig. Topographic surveys were carried 
out determine the elevation river banks and marsh lands. Hydrometric 
surveys were carried out determine the distribution flow during low, 
normal, and flood stages the river system. Sediment sampling was carried 
out various locations the delta determine where sediment was deposit- 
for certain flow conditions. Historical information was collected from 
journals early explorers and from the memory older residents. 

Based the foregoing studies, the following conclusions were reached: 
the present time, about 2,500 acre-ft. per year. sediment load carried 
the Old Channel and remains the river channel well past The Pas. The 
remaining 9,500 acre-ft. carried the New Channel and nearly all de- 
posited the Cumberland Lake Delta. Only 1,000 acre-ft., mostly wash load, 
leaves Cumberland Lake and carried Cedar Lake. Full development 
the delta area includes reclamation all land south the Old Channel and 
Saskatchewan River, and will divert the Old Channel into the New Channel. 
then, the New Channel will carry the full 12,000 acre-ft. sediment load. 
Twenty years later, the bed-material deposits will have reached the outlet 
Cumberland Lake and increasing quantities bed material will begin flow 
down the Saskatchewan River. Thirty years later, the Saskatchewan River, be- 
tween Cumberland House and The Pas, will carry its full capacity bed ma- 
terial load and deposition will begin afterwards, starting Cumberland House 
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and gradually extending downstream direction. The foregoing sedimentation 
graphically illustrated was estimated that 100 yr. from now the 
extreme flood levels the Cumberland House would raised about ft., 
The Pas area 1/2 ft., and near the apex the delta area negligible 
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amount. From then on, the extreme flood levels throughout the whole delta 
area upstream The Pas will continue rise from ft. per 100 yr. 
Temporary and local deviations from this average rise may occur due ab- 
normal climatic conditions due avulsions the river channels. 
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TRAP EFFICIENCY RESERVOIRS, DEBRIS BASINS, AND DEBRIS 


Charles Moore,! ASCE, Walter Wood,2 ASCE 
and Graham 


SYNOPSIS 
attempt has been made review and summarize available infor- 
mation pertaining totrap efficiency and make available engineers 
and the general unified form. additional data are collected 
these and other projects, hoped that will added that pre- 
sented herein that the subject trap efficiency may advanced 


much its importance warrants. 


INTRODUCTION 


Trap efficiency reservoirs subjecton which lackof readily 
available use engineers the design water supply reser- 
voirs, flood control structures, and other kinds impoundment diversion 
structures. Trap efficiency expressed the ratio between sediment ac- 
cumulation and sediment outflow. One objective the Task Force Reser- 
voir Sedimentation, ASCE, has been make summary review the knowledge 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. February, 1960. 
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now accumulated, but not available unified formto the large. 
The authors this paper have attempted this. 

appears that the most authentic information available print paper 
Gunnar Brune.4 This paper covers the normal ponded reservoirs very 
thoroughly this date. summary the paper made herein, including 
the more important findings and results. addition, brief review made 
the cooperative studies trap efficiencies upstream floodwater retarding 
structures currently being made cooperatively the Geological Survey 
and the Soil Conservation the work being done the Los 
Angeles County Flood Control District. There may other published infor- 
mation which the authors are not aware. 


SUMMARY PAPER BRUNE 


pointed out Brune, the trap efficiency reservoir depends upon 
number factors. Among these are the ratio between the storage capacity 
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FIG. 1.—TRAP EFFICIENCY RELATED CAPACITY-WATERSHED RATIO 


and inflow, age the reservoir, shape the reservoir basin, the type out- 
lets and method operation, the size grading the sediment, and the behavior 
the finer sediment fractions under various conditions. number attempts 
have been made tocorrelate trap efficiency with oneor more these factors. 

One the earlier studies trap efficiency was made Brune and 
who developed curve relating percentage eroded soil caught 
reservoir with capacity mile drainage area. The trap-efficiency 
values given this method, however, are necessarily low because they are 


4“Trap Efficiency Reservoirs,” Brune, Transactions, American Geo- 
physical Union, Vol. 34, Number June 

Considerationof Factors Influencing Reservoir Sedimentation the Ohio Valley 
Region,” Brune and Allen, Trans., Amer. Geophys. Union, 22, pp. 
649-655, 1941. 
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based upon rates erosion which, pointed out Glymph, Jr.,6 are 
higher than rates sediment production the reservoir. Glymph found that 
the rate sediment productionis inversely proportional drainage area and 
usually lower than the rate erosion because the deposition some 
the eroded material such colluvium, flood plains and stream channels. 

The ratio between storage capacity and inflow has been expressed gener- 
way the capacity-watershed (C/W) ratio. first developed 
acurve relating C/W trap efficiency. The curv2 shown Fig. 
with some additional records, represented 


where the reservoir trap efficiency percentage, and C/W the reser- 
voir capacity acre feet per square mile drainage area. There con- 
siderable spread points, may noticed. 

The main reason for the spread points the curve that reservoirs 
having the same C/W ratio may have very different capacity-inflow ratios, 
previously stated Gottschalk:8 


watershed ratio, yet not receive enough inflow any one year cause 
water discharged over the spillway. the volume 
mean annual flow from watershedof equalsize humidarea may 
equivalent times that reservoir having the same capacity- 
watershed ratio. the dried region 100% the incoming sediment load 
trapped, whereas humid area possibly only 70% trapped.” 


Allen first introduced the storage, capacity-inflow, ratio 1914. 
used for determining reservoir requirements, however, and not in- 
dex sediment trap efficiency. 

For individual reservoirs, curves can drawn correlating trap efficiency 
with detention time days. Whitney Borland has such curve 
for Imperial Dam Reservoir, shown Fig. Such curves are quite satis- 
factory for specific reservoirs, since other factors such sediment charac- 
teristics, shape method operation tend remain constant. 
Detention time may defined the period that the inflow remains the 
reservoir before being released downstream. 

took into account both detention time and velocity flow 
through the reservoir. hasdeveloped “sedimentation index” which repre- 
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sents the period detention divided mean velocity. His curve, relating 
the sedimentation index for several TVA reservoirs, shown 
Brune’s paper. While this curve very satisfactory where the data are 
available, such information period detention and mean velocity not 
readily available for most reservoirs. 

Brune’s study thorough search was made for all reliable records 
reservoir trap efficiency, canvassing all agencies the Federal Inter-Agency 
River Basin Committee known interested the subject. Information was 
also gathered capacity, annual inflow, shape reservoir basin, type out- 
lets and method operation, observed gravity underflows density” cur- 
rents, and any other pertinent data. 

Some forty-four records were gathered and are summarized Table in- 
cluded herein. these, forty are for normal ponded reservoirs, with infor- 
mation the effect sluicing and venting operations three. Two records 
are for desilting basins, and two for semi-dry reservoirs. 


TRAP EFFICIENCY 


From 'cansactions, American Geophysical Union, Yol. 3%, Ho. 3, 1953. 


DETENTION TIME IN DAYS 


FIG, 2.—TRAP EFFICIENCY RELATED DETENTION TIME, IMPERIAL 
DAM RESERVOIR, YUMA, ARIZONA 


pointedout Brune, possible study trapefficiency reservoirs 
number different methods, and found necessary use nearly all 
them. Those interested the details are referred Brune’s paper. 

taking suspended-load measurements downstream from reservoir 
necessary near the dam possible, because downstream from 
reservoirs streams tend pick new sediment load quickly scouring 
and degrading their channels. For example, Brown has pointed out,12 
the average sediment concentration the Colorado River below Hoover Dam 
1947 increases from 0.006% 12.8 miles below the dam 0.091% 98.6 
miles below the dam, times. Such change can cause grave errors 
computations trap efficiency suspended-load measurements far down- 
stream from the dam are used. 

computing trap efficiency also necessary take dredging into ac- 
count. Since dredged areas reservoir usually fill with sediment again 


“Sediment Transportation,” Brown, Chapter XII Hydrau- 
lics,” Hunter Rouse, John Wiley and Sons, New York, Y., Aug. 1950b. 
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short time, dredging tends increase the amount sediment trapped 
reservoir. 

Table reservoirs and ponds allsizes were used, ranging froma farm 
pond with 0.038 mile drainage area Imperial Dam Reservoir, draining 
184,600 miles. The chief criteria selecting reliable records were accu- 
racy with which trap efficiency could determined and the availability in- 
forination the various factors which affect trap efficiency. Gottschalk 
stated,8 the laws sediment deposition are the same for all types reser- 
voirs, including stock ponds. Likewise, the trap efficiency reservoirs af- 
fected the same factors, regardless the size the reservoir. 

has long been that capacity-inflow (C/I) ratio would 
more accurate indexof trap efficiency than the capacity-watershed (C/W) ratio 
which has heretofore been widely used. The C/W ratio must used within 
definite hydrologic regions, and notas indexof the country 
whole, for, obvious, with the same C/W ratio the trap efficiency will 
increase the runoff per unit area decreases. 

Table Brune’s paper reproduced herein and the records have been 
plotted shown Fig. Brune found that correlation for normal ponded 
reservoirs (conventional reservoirs distinguished from desilting basins and 
dry reservoirs) operated without any special efforts sluicing sediment, 
much better using the C/I ratio than using the C/W ratio. 

Records 32-37ain Table Fig. represent reservoirs suchas Fort 
Peck Reservoir and Lake Mead which have C/I ratios greater than1.0. 
other words, these reservoirs have capacity larger than the annual runoff 
from their watersheds. When the effects evaporation and seepage are also 
considered, obvious that water will rarely discharged over the spill- 
ways such reservoirs, and that trap efficiency must therefore close 
100%. 

The C/I ratio thus provides means differentiating between the “hold- 
over storage reservoirs” and “seasonal storage reservoirs” described 
Thomas Maddock, Jr.13 Reservoirs with C/I ratio 1.0 less may 
classed seasonal storage reservoirs, and those with ratio greater than 
1.0 holdover storage reservoirs. 

Reservoirs very different C/W ratios and annual inflows, yet have 
the same C/I ratio and the same trap efficiency, other factors being equal. 
For instance, reservoir the humid part the country may have C/W 
ratio 100, annual inflow 1,000 acre per sq-mile, and C/I ratio 
0.1. Another reservoir arid region may have C/W ratio 10, an- 
nual inflow 100 acre per mile, and the same C/I ratio Both 
reservoirs, assuming they are the normal ponded type, will have roughly the 
same trap efficiency (around 87%, from Fig. 3). 

Theoretically, reservoir can have trap efficiency zero only when its 
capacity reduced zero the runoff approaches infinite quantity. Also, 
trap efficiency 100% theoretically possible only when the reservoir can 
contain the inflow long enough allow all sediment out. inspection 
Table and Fig. shows that there instances record trap 
efficiencies zero and 100%. 


“Reservoir Problems with Respect Sedimentation,” Thomas Maddock, Jr., 
pp. 9-14, Proc., Federal Inter-Agency Sedimentation Conference, Washington, C., 
1948. 
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Williams Reservoir, power reservoir near Williams, Inc., had trap ef- 
ficiency less than zero during the period measurement between 1930 and 
1939, thus revealing that scouring actually had occurred. Three instances 
trap efficiency equal 100% were measured Lage Pond, Fort Peck Reser- 
voir, and Denison Reservoir. Details these records may foundin Brune’s 
paper. 

Table and Fig. show the effects sluicing and venting. Arrowrock 
Reservoir, Idaho, sluicing estimated have reduced the trap efficiency 
from 93.0% 90.3%, and have increased the sediment outflow from the reser- 
voir 39%. Likewise, Conchas Dam, M., sluicing operations decreased 
the estimated trap efficiency from 97.3% 95.8%, and increased the sediment 
outflow from the reservoir 56%. Those who are interested the details 
these and other sluicing and venting operations are referred the original 
paper. 

Brune, appears that sluicing and venting operations have 
much better chances accomplishing their purpose when they are timed 
intercept gravity underflows they reach the dam. gravity underflows are 
not present, such operations should least timed meet the higher sedi- 
ment concentrations brought flood flows. likely that accurate timing 
intercept gravity underflows can quadruple the amount sediment 
lost from the reservoir. 

All the foregoing pertains normal ponded reservoirs. word now 
order about desilting basins and dry reservoirs studied Brune. The 
trap efficiency desilting basins governed largely their shape. Hadley 
Creek New Desilting Basin, and Fig. covers 280 acres but has 
average effective depth less than ft. normal ponded reservoir with 
the same C/I ratio would have trapefficiency about 85%. This basin, how- 
ever, traps 98.8% the sediment reaching it, largely because its shallow- 
ness and the small interval time required for the sediment reach the 
bottom. Other factors, such flocculation, the coarse nature the sediment, 
and the location the basin point where the stream gradient decreases 
abruptly, have their effect increasing the trap efficiency. 

even more efficient type desilting basin exemplified the All- 
American Canal Desilting Basins, Ariz. The method operation, addition 
the factors previously mentioned, causes these basins have high trap 
efficiency. normal ponded reservoir with C/I ratio low these basins 
could expected very low trap efficiency, yet the basins trap 91.7% 
the sediment reaching them. 

Thus, according Brune, desilting basins may have trapefficiencies above 
90% much lower range C/I ratios than for normal ponded reservoirs. 

Although information available efficiency truly dry reser- 
voirs, two reservoirs which have been operated semi-dry condition pro- 
vide indicative data. 

The Senecaville Reservoir, Ohio, flood-control reservoir built the Corps 
Engineers, and the John Martin Reservoir, Colo., are reservoirs 
being operated semi-dry condition. The studies the Senecaville Reser- 
voir revealed the trap efficiency 48.4% when operated with the gates wide 
open and 94.3% with the gates closed. The trap efficiency the John Martin 
Reservoir was 62.2%, unusually low figure for reservoir with such high 
C/I ratio (1.71). 

Thus apparent that even holdover reservoirs like John Martin, oper- 
ated allow large flows water pass through the dam, may have 
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TRAP EFFICIENCY 


trap efficiencies the 60% range rather than about 90% would expected 
with normal operation. Truly dry reservoirs, such the Miami Conservancy 
District reservoirs, undoubtedly have much lower trap efficiencies, probably 
the 10% 40% range, depending chiefly upon the C/I ratio. 

The more important conclusions reached the study Brune have been 
summarized the foregoing portions the paper. However, those who are 
interested more details are referred the original paper. There are nu- 
merous references listed the end his paper. 


UPSTREAM FLOODWATER RETARDING STRUCTURES 


pointed out previously, the only known published information available 
trap efficiency that Brune, which has been summarized. Since much 
this information applied storage reservoirs, not known whether 
applies also the type floodwater retarding structures (particularly with 
dry pools) being constructed the Soil Conservation Service inthe watershed 
protection program. For this and other reasons, the trap efficiency studies 
obtain additional information for design were initiated the Soil Conservation 
Service under authorization Public Law 566 cooperation with the 
Geological Survey. The studies are being executed the Geological Survey 
with funds transferred the Soil Conservation Service. 

The effectiveness floodwater retarding structures trapping sediment 
from storm runoff has been observed only general way. The objective 
the trap efficiency studies provide planning data for the design such 
structures the range watershed conditions exhibited the United States, 
delineate the principal factors that influence trap efficiency, and study 
the sediment yield and characteristics small watersheds. The selection 
individual projects has been made jointly representatives the Geological 
Survey and the Soil Conservation Service and has been guided the relative 
need for design data, the availability supplementary hydrologic data and 
the diversity watershed characteristics. 

Before proceeding further with the discussion the cooperative trap- 
efficiency studies might well discuss briefly the nature and general 
design upstream floodwater retarding structure. Fig. section 
typical structure. The structures are rolled earth embankments with earthen 
vegetated spillways. They are designed without conservation pools. Their de- 
sign does provide storage space for sediment, the needed space being deter- 
mined sediment source survey the watershed each structure. the 
design for this storage, sediment rates based present conditions the water- 
shed for ten years and conditions for are used where the 
structure designed operate full effectiveness for 50-yr period. De- 
pending the character incoming sediment, from 10% 30% the sedi- 
ment storage allocated the flood pool, the remainder being provided be- 
low the inlet the principal spillway the sediment pool. 

The principal spillway floodwater retarding structure usually consists 
reinforcedconcrete drop inlet with heavily reinforced concrete pressure 
pipe, having minimum diameter in., laid concrete cradle through 
the embankment. 

The principal spillway automatically releases water from temporary stor- 
age the detention pool controlled release rate. Its discharge capacity 
determined from consideration the various physical and economic factors 
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associated with the planning single dam system retarding dams. 
Variations discharge capacity influence detention storage time, and corre- 
latively the trap efficiency. 

sluice gate may constructed the inlet the principal spillway 
permit the release water stored below the crest pending the accumulation 
sediment. This essential some areas satisfy downstream water 
rights. 

The capacity the detention pool such, that flow through the emergency 
spillway occurs not more often than once for any structure. Many 
such structures are designed with sufficient detention storage accommodate 
flood 100-yr more estimated frequency before the emergency spillway 
brought into operation. 

The sediment and flood pool volumes vary with the size drainage area, 
runoff characteristics, sediment yield and other factors. general, however, 
the sediment pool storage lessthan 25% the total available storage behind 
the structure. the Southwest area the average volume sediment space 
209 acre compared 1,069 acre total storage, average for 146 
structures. The character sediment, character native water, and detention 


storage time appear the most significant parameters influencing trap ef- 
ficiency. 


FIG, 4.—TYPICAL SECTION FLOODWATER RETARDING STRUCTURE 


much for the general design the floodwater retarding structures 
might affect trap efficiency. 

The operating procedures are adjusted the particular characteristics 
each structure. Twelve structures eleven states were selected repre- 
sentative for trap-efficiency studies.14 Information the operation the 
structures June 30, 1956 covers only the first ten structures. The last 
two structures were added after that date. 

The station name, location, date established, area watershed, capacity 
the detention pool (from sediment pool elevation the earth spillway ele- 
vation), and the detention time each project are shown the report the 
order their installation. The procedures each project location are de- 
signed provide: (1) qualitative data the character water and sediment 
flowing into the reservoir; (2) quantitative data the amount and properties 
sediment retained the pool area; and (3) quantitative data the amount, 
rate, and properties sediment discharged from the structure. Occasional 
samples inflowing storm runoff, the order about per collect- 


Report, Investigation Trap Efficiency Floodwater Retarding Struc- 
tures, July 1955 June 1956,” Geological Survey Soil Conservation Service, 
Feb. 11, 1957. 
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describe the concentration and size distribution inflowing sediment and 
pertinent chemical properties the inflowing water. The Soil Conservation 
Service has each instance the responsibility making periodic pool sedi- 
mentation surveys determine the amount, distribution and characteristics 
the deposits. The frequency pool surveys determined the rate 
accumulation and the occurrence unusual flood events. The amount sedi- 
ment discharged from the structure measured sampling frequency 
sufficient todelineate continuous recordof sediment concentration with time. 
Sampling frequency ranges from hourly intervals during critical periods 
storm discharge weekly intervals during periods base flow discharge. 
The size distribution and flocculation characteristics the sediment, and perti- 
nent chemical properties the water, are determined about ten times per 
year. 

Supplemental data relating precipitation, runoff, inflow-outflow relation- 
ships, physical and cultural characteristics the watershed and physical 
characteristics the structure are available each project watershed from 
concurrent investigations the Soil Conservation Service and Geological 
Survey. 

Operation five the projects was begun the interval between January 
and June 1956. Study the two reservoirs Nebraska unique that struc- 
ture No. within the watershed structure No.1. the present time only 
brief information covering the results the studies available. The estab- 
lishment and testing relationships used computing water and sediment 
discharge customarily cause more lag between observation and computation 
data for the initial period than for subsequent periods operation. 

Since results measurements June 30, 1956, the latest available in- 
formation some the structures, are meager due the shortness the 
the few runoffevents since initiation the studies, they are not be- 
ing included. 

However, the numerous storm events that occurred during March, April and 
May 1957 Texas, Oklahoma and Arkansas provided considerable information 
and data relative inflow and outflow for the floodwater retarding 
structures which trap-efficiency studies are being These data are 
being analyzed and tabulated the Geological Survey. 

Although sufficient time has not elapsed and the storm events above the 
structures have not been numerous enough give suitable data which tode- 
sign, appears that goodstart has been obtain reliable and accurate 
data needed for evaluation trap efficiency for structures this type. Until 
such time the data does become available the Soil Conservation Service, 
among others, will necessity have adapt the presently available information 
for other types reservoirs the type structure being considered. this 
respect, indications are that, for the floodwater retarding structures the 
type being designed the Soil Conservation Service and having sediment pools 
such those constructed the Southwest area, 85% 95% the incoming 
sediment trapped. This appears reasonably conservative estimate, 
but will revised more accurate data become available. 


“Sedimentation Studies the Soil Conservation Service the Western Gulf 
States,” Renfro and Moore, Paper presented National ASCE Conven- 
tion, Hydraulics Division, Jackson, Miss., Feb. 20, 1957, and published the Journal 
the Hydraulics Division, ASCE, Oct. 1958, Paper No. 1806. 
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RESERVOIRS, DEBRIS BASINS AND DEBRIS DAMS 


Before reviewing and summarizing the work being done the Los Angeles 
County Flood Control District the control debris-laden streams, might 
helpful know something the background the debris basin program 
the Los Angeles area. The following information taken from letter report 
prepared December 1952 member the District’s staff the subject 
“Control Debris Laden Streams.” 

Southern California has experienced severalsevere floods the past which 
have produced heavy from the mountains. One notable example 
the flood resulting from the storm December 31, 1933 January 1934. 

brush fire November 1933 burned over area 4,500 acres the 
mountains north Crescenta and Montrose about miles north Los 
Angeles. The mountain drainage this area has very rugged topography 
rising some 2,500 horizontal distance 1-3/4 miles. The average 
slope the debris cones the foot the mountains between and 10%. 
Individual watersheds these mountains range from 0.2 1.6 miles 
area. 

During the storm mentioned, more than in. rain fell. This amount 
rain falling steep, recently burned-over areas produced great quantities 
produced from the 4,500 acres during this storm. Although velocities mud 
flows were sufficient move many large boulders, tons more weight, 
the bulk the debris produced was probably less than 1/4 in. size. 

While heavy rainfall recently denuded surface will produce great 
quantities debris,a major storm canproduce similar quantities un- 
burned saturated area. For example, the storm March 1938 fell fairly 
well saturated ground and produced about 72,000 debris per mile 
from the same area. One area 0.27 miles produced debris rate 
120,000 per mile. 

the time the flood January 1934 gravel pit existed near the 
mouth Haines Canyon. This pit had been provided with temporary spill- 
way investigate its action debris basin,collecting the debris and al- 
lowing desilted water pass down the channel. The debris damage down- 
stream from this notas severe the similarly protected. 

With this example before it, the Los Angeles County Flood Control District 
1934 embarked program debris basin construction. Later the 
Army Engineers undertook the construction additional basins the Montrose- 
Crescenta area. 

Essentially, this system basins located the mouths potenti- 
ally hazardous canyons andof concrete-lined lateral channels leading from the 
basins major collector channel. The lateral channels being located 
debris cones, are naturally steep grades and consequently carry the flow 
high (super-critical) velocities. Careful attention must therefore given 
super-elevation curves and oscillation the water prism thereby creat- 
ed. Rectangular reinforced concrete channels have been found the most 
satisfactory for these conditions flow. 

mentioned, the New Year’s storm January 1934, produced about 
87,000 debris per mile watershed. This led the adoption 
design capacity value 100,000 per mile for debris basins the 
Montrose-La Crescenta area. This has been the accepted design value when 
permanent installation considered. 
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However, since debris production dependent upon number variables 
such rainfall intensity, slope watershed, vegetation and soil conditions, 
has been the experience the District that the rate debris production per 
square mile for areas other than the Montrose-La Crescenta area can range 
between 40,000 and 10,000 yds per mile. view this wide variation, 
some means estimating previous amounts debris flows from given 
watershed provides the most satisfactory means establishing production 
rates. present there completely reliable means using known de- 
bris for given watershed determining the production rates 
for different watershed, except general comparison topographic fea- 
tures and rainfall characteristics. 

The basins themselves have been two types, (1) excavated, and (2) non- 

The excavated obtained digging pitof the required capacity and 
using the material removed form dam suitable curvature around the 
downstream the basin. this dam adequate concrete spillway struc- 
ture height constructed, generally line with and having its 
foundation the existing streambed. Several uncontrolled discharge ports 
are generally incorporated the spillway allow for some release during 
storms, and partial post-storm drainage. essential that intake struc- 
ture provided from the normal streambed bottom the excavated pit. 
Otherwise, channelerosion will occur the inlet tothe basinand progress up- 
stream. line connecting the lips the inlet and spillway structures gener- 
ally slope that about one-half two thirds the natural slope the 
stream channel. These basins are particularly appropriate where complete 
desilting, least the ports, desirable. 

The unexcavated debris basinis similar small storage reservoir, with 
few variations. The dam generally compacted earthfill with heavy rock 
armoror reinforced concrete its upstream slope. For purposes estimat- 
ing the storage capacities, the spillway lip elevation that will create 
the desired debris storage capacity. The debris slope assumed about 
one-half two-thirds the natural channel slope upstream from the spillway 
lip. inlet structure required for unexcavated basin. provision 
generally made allow for some release during storms. This may done 
means ports through the spillway some other outlet structure. These 
basins lend themselves partial desilting, that is, the retention debris 
such size only would detrimental downstream 
from the basins allowed discharged therefrom. This may accom- 
plished designing the width and slope the basin soas cause the desired 
segregation debris take place automatically. 

has been the policy the Los Angeles Flood Control District limit 
debris accumulation 25% the original basin capacity before debris re- 
moval considered necessary. There fixed rule for removal, however, 
for rather flexible basis has been used the operation and maintenance 
the debris control system. The allowable debris accumulation considered 
the the probable debris production and the condition the watershed 
relative vegetation and geological characteristics. 

Debris disposal presents major problem rapidly growing southern 
California. Generally speaking, the accumulated debris has agricultural 
value spite the very high percentage fine material. After number 
years operation debris control system with the consequent debris re- 
moval, the nearby disposal areas are exhausted and more remote sites must 
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found. This makes mechanical debris removal costly procedure. would 
desirable segregate the bulkof the fine debris and cause discharged 
with the water, previously discussed. This type debris removal elimi- 
nates the the accumulated debris and minimizes interference 
with the normal replenishment sand for the beaches southern California. 

general, floating debris does not present great problem for the debris 
basins. The watersheds are usually quite small and the volume this type 
debris not large. However, until the water surface the basin rises above 
the ports, there always the possibility clogging. Patrolmen have been 
able keep the ports clear more less satisfactorily during storms. The 
District has been experimenting with inclined made from railroad 
rails and spaced roughly one-half the minimum dimension the ports. The 
been satisfactory far, though incomplete. When the water over- 
flows the spillway lip the floating debris allowed pass down the chan- 
nel. the end the runoff season when the basin has been allowed drain, 
the remainder the floating debris accumulates the sides and bottom. 
then either collected and burned removed from the basin trucking. 

times itis necessary construct temporary debris basins and channels 
before sufficient funds are available build permanent structures. Under 
these conditions the design capacity, 100,000 per mile, sometimes 
reduced per mile, less. Such basin might have asphalt- 
lined wooden inlet and spillway structures, well asphalt-lined dis- 
charge channel. The use asphalt replace concrete requires caution, par- 
ticularly where impact due coarse debris and boulders concerned. Asphalt 
will not stand under such punishment. will, however, resist abrasion due 
ordinary drag bed load debris, the maximum size which the 
order three inches. 

Necessarily, any flood protection must based definite plan con- 
trol. Portions area might require storage reservoirs for debris control 
and flood flow regulation while smaller watersheds may needonly system 
debris basins for satisfactory protection. 

Since December 1952 the most important modification the early debris 
basin concept and design has been the elimination dead storage and the 
open port outlets which were provided addition the spillway many 
the debris basins, and the substitution outlet tower. This located 
the low point the bottom the debris basin offset from the line flow be- 
tween the inlet and the spillway and serves the entrance conduit which 
passes under the dam and discharges the channel downstream. This outlet 
design accomplishes two things: (1) permits storm flows which are rela- 
tively free debris pass through the basin with minimum ponding, 
thus leaving basin capacity for debris storage, and (2) serves trash 
rack. These towers have been constructed circular form 5-ft inside di- 
ameter reinforced concrete reinforced gunite. They have 6-in. wall 
thickness with series round perforations four inches diameter the 
outside and five inches diameter the inside face the tower. The tower 
sufficient height project foot above theoretical debris slope, 
equal 75% the original streambed slope. The outlet pipe normally 
36-in. reinforced concrete pipe with slope not and capacity 
least 150 cfs for the maximum expected head water the basin. 
means controlling orifice the base the tower, flowunder pressure 
the outlet prevented. Other modifications design have specific bearing 
the subject trap efficiency and need not discussed here. 
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presenting information relative the reservoirs, debris basins and de- 
bris dams the Los Angeles County Flood Control District, appeared that 
perhaps the best way would formatof Table Brune’s paper.4 
This procedure was followed insofar the data permitted. Tables entitled 
“Reservoir Debris Production Data,” Table and “Debris Dams and Debris 
Basins—Debris Production Data,” Table are included herein. 

From the data available was not possible include three columns which 
were included Brune’s paper. These were (1) average annual inflow, (2) C/Ib, 
efficiency. This omission may explained stating that relative 
debris basins and the majority the debris dams, continuous inflow and 
outflow records are not available. fact, very few records inflow and out- 
flow are available. Due the lack inflow records, the values for and 
trap efficiency could not computed. However, believed that nearly 
all cases the trap efficiency, least course-grained debris, approaches 
100%. With the exception five the basins, namely, Aliso Creek, Scholl, 
Snover, Ward Canyon, and Verdugo, the material trapped, will later ex- 
plained. the abovementioned basins, the silt has been allowed accumulate 
the level outlets. Much the fine material and that which suspen- 
sion has been allowed through. This has been somewhat the nature 
experiment, but far appears successful plan operation for 
these basins, least during seasons where the watersheds are not burned off 
and subsequent floods occur. the case Scholl and Ward Canyon, flash- 
boards are placed the outlet notches regulate some extent the level 
debris the basin. the other basins, except the case sizable floods 
which fill the basins with sediment, all material and most the inflow 
trapped during given storm. If, subsequent storm, found that less 
than 75% storage capacity remains, the material removed shovel and 
trucks nearby disposal area. 

The omission the three-above-mentioned columns the tabulation for 
flood control dams due the fact that although the District does have in- 
flow and outflow records, sediment generally goes through the outlets only 
during periods controlled sluicing. Under certain other conditions such 
during flood releases, small quantities silt and sand may pass 
through the valves. Normally flood releases are made only whenan appreciable 
head exists the valve. Other releases are made subsequent storm peri- 
ods and therefore after most the suspended material has settled out. 
stated above, the District believes that trap efficiency for the reservoirs 
generally approaches 100%, even during periods when both inflow and outflow 
take place. the reservoirs dual purpose flood controland water 
conservation, water stored limits compatible with flood operation. 

some the reservoirs where conditions permit, capacity regained 
the hauling some the silt accumulation permittees. general, this 
does not constitute large storage reclamation gained sluicing, 
may noted from the table. Sluicing, course, dependent upon sufficient 
residual flows. This often problem, the reservoirs may not drawn 
down due conservation commitments until late fall when inflows 
are low. 

study has recently been completed the District’s Foundation and Test- 
ing Division, relative the quality debris excavated fromthe debris basins 
and removed disposal areas; the latter have been provided and set aside 
the District exclusively for disposal debris. Inasmuch the size ma- 
terial trapped basins interest and has bearingon the sub- 
jectof trap study was made grading analyses for material from 
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q LOS ANGELES COUNTY FLOOD CONTROL DISTRICT 
FOUNDATION AND TESTING DIVISION 
Hole No. Composite of 5S Pits Lab. Tested by P.P.ac.L Date 5/57 
q Sample Wo. Depth Percent Rock 45.3 Percent Sand 48.4 
Loc. Disposal Area Elev. FL OM. Rock 7.70 Sand 2.17 
ASTM SIEVE NUMBER 
100 200 100 50 30 16 8 3/8° 3/8" 
iy 
Comparison of Debris Material with 
80 Pit Material from San Fernando Valley 
—— Rock Producing Companies am 
= 
= = 
= 
Livingston 
Consolidated 
CLAY AND 
(LIMITS PER A.S.T.m. C 125) 
FIG, 
LOS ANGELES COUNTY FLOOD CONTROL DISTRICT 
FOUNDATION AND TESTING DIVISION 
Project Halis Debris Basin GRADING ANALYSIS Sampled by o.2, Date s/s? 
Hole Mo. Composite of 3 Borings Lab. Tested by P.P.&C.L. Date 5/57 
Sample Mo. Depth Percent Rock 24.0 Percent Sand 64.4 
Loc. Disposal Area Elev. F. Rock 8.20 F. M. Sand 2.36 
ASTM SIEVE NUMBER 
100 200 100 50 30 16 8 4 3/8" 3/4" 13" 3° 
Comparison of Debris Material with : 
90 Pit Material from San Fernando Valley 
—— Rock Producing Companies 
= 
= 
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Livingston 
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basins including Cooks, Dunsmuir, Halls, Pickens, Ward, and West Ravine. 
The grading analyses for the various basins were compared with pit material 
from the San Fernando Valley rock producing companies. Figs. are 
typical gradation curves for material from the Dunsmuir and Halls debris 
basins. 

The rock content the debris material found the various disposal sites 
varied from 10% 45%, with average 29%. This considerably less 
than the average about 58% rock found the commercial material. 

general, the gradation the rock was considerably finer than that the 
commercial sources. Excluding the material from West Ravine, Halls, and 
Aliso-Wilbur debris basins, the debris rock compared favorably with that 
the commercial sources relative abrasive loss. 

The debris material from the various disposal sites contained between 
48% and 74% sand, with average 63%. This greater than the average 
40% sand found the commercial material. The quality the debris sand 
compared favorably with the commercial material. 

order determine more accurately the degree trap efficiency the 
Los Angeles County Flood Control District structures, will necessary 
obtain more data directly related the subject. These data should include 
continuous inflow and outflow records and measurements sediment coming 
into and passing through the structures. such data were available, the trap 
efficiency could computed C/I ratio basis and perhaps correlation 
could made for these structures has been done for normal ponded reser- 
voirs. 
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STRATEGIC ASPECTS URBAN FLOOD PLAIN 


Gilbert 


SYNOPSIS 
The paradox continuing flood control and rising flood losses ex- 
plained part elements entering into decisions change urban 
flood-plain occupance, 1936-1957. Curbing the mounting toll calls for 
public action broaden the range choice open managers flood- 


plain properties adjusting flood hazard. 


INTRODUCTION 


striking feature national efforts flood control the United States 
that the more than two decades since national policy was launched 1936, 
the mean annual toll flood losses has continued rise. While least bil- 
lion dollars have been spent for engineering works reduce and control 
floods, the economic losses from occasional flood disaster have mounted. 
understand this paradox necessary examine what happening 
the human occupance flood plains. find out the full import changes 
flood-plain occupance essential identify those factors which enter into 
both public and private decisions flood-plain use. 

This paper outlines briefly the broad problem flood-loss reduction the 
United States, findings from recent study changes urban oc- 
cupance flood plains, and then suggests three strategic aspects such oc- 


cupance that seem affect the results any further engineering efforts 
flood control. 


Note.—Discussion open until July 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. February, 1960. 

Presented the July, 1959 ASCE Hydr. Conf. Fort Collins, Col. 

Dept. Geography, The Univ. Chicago, Chicago, 


a 


The broad problem flood-loss reduction that the rate which flood 
losses are being eliminated construction engineering land-treatment 
works about the same magnitude the rate which new property be- 
ing subjected damage. Even though the data which estimates national 
flood loss are based inaccurate precise comparison 
the two rates, clear, both from the aggregated damage statistics and from 
the record selected flood plains, that the heavy investments flood protec- 


THE NATIONWIDE FLOOO 


PROBLEM 


RESIDUAL AFTER 

AUTHORIZED PROGRAM 

COMPLETE 


se 
WORKS 


UNDER CONSTRUCTION 


EXCLUSIVE DAMAGES SMALL UPSTREAM 
TRIBUTARIES ESTIMATED DEPT AGRICULTURE 
MILLION ANNUALLY. 


FIG. 1.—GRAPH PRESENTED THE CORPS ENGINEERS THE HOUSE 
COMMITTEE PUBLIC WORKS 1957 SHOWING ESTIMATE 
MEAN ANNUAL FLOOD-DAMAGE POTENTIAL AND DAMAGES 
PREVENTED. (FROM 85TH CONGRESS, 1ST SESSION, HOUSE 
COMMITTEE PUBLIC WORKS, PRINT NO. 1957). 


tion have effectively curbed the losses many areas but that new damage po- 
tention being built the same time. 

Studies seventeen selected urban areas having flood problems reveal 
general and persistent encroachment urban structures upon the flood plains 
during the period 1936-1957, even areas which there was net decrease 
total population. They show certain distinctive patterns encroachment, and 
indicate highway construction and flood-control works astwo major stimulants 
growth. The situations Boulder, Colo. and Denver, Colo., illustrate some 
the major findings. 
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the light that experience argued that least three aspectsor ur- 
ban occupance not ordinarily considered flood-control plans the past must 
taken into account the future the tide rising flood losses 
turned. First, must recognized that engineering works are only one the 
possible human adjustments flood hazard. Second, must seen that the 
complex elements entering into decisions future occupance flood 
plains includes many considerations addition the traditional one cost- 
benefit evaluation. Third, the range choice now permitted property man- 
agers dealing with the flood hazard restricted that radical changes 
must made public policy broaden the range choice among possible 
adjustments and assure full appraisal each choice. 


TRENDS FLOOD LOSSES 


There are two major sources for national estimates flood losses. The 
Corps Engineers issues, from time time, estimate the total poten- 
tial flood losses and the losses that have been prevented will pre- 
vented works constructed planned the Corps. sample graphical 
statement such estimate for 1954 given Fig. that year the 
Corps found potential mean annual loss 964 million dollars.2 (All esti- 
mates are adjusted 1957 price levels.) After all works then authorized had 
been completed there would remain approximately 210 million dol- 
lars that would not prevented. These losses are computed separately for 
each the major hydrographic areas the country (Fig. 2). 

The Weather Bureau issues the other series flood-loss statistics, having 
published each year since 1903, estimate reported losses for each its 
districts. The national totals for 1902-1955 are shown the lower part the 
bars the graph Fig. These totals have been adjusted 1957 price 
levels, and the full height the bar shows the totals comparable basis. 

attempt will made here appraise the validity and the discrepancies 
among these two sets estimates and other less comprehensive series. This 
has been done separate Aggregate data alone are unsatis- 
factory and need checked more detailed fashion. sufficient here 
point out one characteristic which the two national estimates have com- 
mon. pronounced upward trend the size annual 
since 1936, when the first national flood control legislation was enacted. The 
estimated mean annual losses 1936, were about 212 million dollars. The 
Corps estimate remaining annual losses after allowing for works completed 
was 444 million million 1959, (at 1959 prices).4 The Weath- 


Bureau estimates mean annual losses for 1924-1953 were 25% lower than 
those for 1944-1955. 


“House Committee Print,” Congress, House Committee Pub. Works, 1957, 
No. 85th Congress, Sess. 

Urban Occupance Flood Plains the United States,” Gilbert 
White, Wesley Calef, James Hudson, Harold Mayer, John Sheaffer, and 
Donald Volk, Univ. Chicago Geography Research Paper, 1958, No. 57, pp. 1-19. 

“Hearings Civil Works Appropriations Bill,” Congress, House 
tions Committee, 1959, Pt. 18. 
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REASONS FOR MOUNTING FLOOD LOSSES 


One partial explanation for the apparently rising toll recorded losses 
that the inflation dollar values has exaggerated the size recent estimates. 
shown Fig. even when values are adjusted common year there 
still remains clear and sizeable increase. 

second reason advanced for the apparently mounting losses improve- 
ment enumeration and estimating methods. This unquestionably has played 
role Weather Bureau procedures have been made more uniform and 
Corps Engineers studies, have reached areas not previously covered. Some 
who have been studying this problem believe that changes enumeration meth- 
ods may account for much 10% 15% the increase. 

Perhaps greater importance the difference number large 
frequent floods. William Hoyt, F., ASCE and Walter Langbein, F., ASCE 
have shown 35% increase_in their “flood index” from the first yr. the 
century the second There was remarkable bunching rare flows 
the lower Missouri and Northeastern basins during the 1950’s. This hydro- 
logic record may account for much 25% the increase losses. 

Probably the most important reason for the rising trend flood losses 
found the continuing encroachment human occupance upon flood 
plains. This takes the form new structures, changes the intensity 
existing structures, and structures which reduce the hydraulic efficiency 
valley sections increase the hazard affected reaches the stream. 
Although there have been numerous studies flood-control projects there has 
been comprehensive investigation the changes that have been occurring 
the areas subject flood. One recent study conducted University Chicago 
geographers dealt with urban occupance because has high and increasing 
proportion flood losses, compact and more susceptible change. 
shows that there are least 1,020 urban places with populationof more than 
1,000 that have well defined flood problems. That study also gives precise 
picture the changes seventeen urban areas selected for their diversity 
valley section, flood frequency and height, size, population growth, and types 
land use. few findings from those field examinations will summarized, 
and then illustrated from the Boulder and Denver situations. 


GROWTH PATTERNS 


The most evident and widespread trend that observed the urban 
flood plains one growth. every place studied, including several 
which the total population declined during the period, the number 
structural units the flood plain increased.6 The growth rate ranged from 
less than stable city like Wheeling, Va., more than 600% 
rapidly expanding place such Dallas, Tex. 


“Floods,” William Hoyt, and Walter Langbein, Princeton Univ. Press, 1955. 
Changes occupance were measured “structural unit was defined for 
different occupance classes as: Residential, A-a single double-family dwelling, B-a 
multi-family dwelling for 3-6 families, and C-a multi-family dwelling for more than 
families; Commercial, Industrial and Transport, separate structures multiples 


10,000 sq. ft, and open working space multiples 25,000 sq. ft; Public, each separate 
building. 
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most places the commercial and industrial structures grew more 
rapid rate than residential structures, although few the places such 
sectors Los Angeles, Calif., the residential expansion was substantial. 
Public structures most places showed substantial growth rates. 

Certain areal patterns change seem repeat themselves urban flood 

4). One typical pattern found where residential areas bor- 
dering flood plain enlarge moving down nearer the channel. second, 

that which industries already established flood-hazard zone expand 

laterally along the stream bank and right angles away third oc- 
curs where central business district spreads out into the flood plain, follow- 
ing the major traffic arteries. 

There readily discernable relation between rate population growth 
and rate change structures the flood plain: Too many other factors 


New 


gy Qiadestriet Nishwey Reilwey 


COMMERCIAL 


RESIDENTIAL INDUSTRIAL 


FIG. 4.—SCHEMATIC MAPS THREE STRATEGIC LAND-USE ASSOCIATIONS, 
SHOWING PATTERNS GROWTH THE FLOOD PLAIN. 


such urban function and availability landplay part. Itis clear, however, 

that two forms public action have had powerful stimulating effect upon in- 

vasion flood plains. These are the highway program and the flood-control 

program. New highway construction urban areas tends follow the low 

gradients stream valleys and the less densely settled sectors some flood 

plains. has caused substantial removal low quality residential structures 

but also has been inducement for commercial and industrial establish- 
ments follow along into the flood hazard zones. 

might expected, the construction new flood-protection works fre- 

quently has been the signal for accelerated movement into the flood plain. 

Thus, the completion reliable works along the Trinity River Dallas, Tex. 

assured development large-scale commercial district behind the levees. 


is 


But there are certain implications flood-control works that are not 

ous. The lack protection not necessarily deterrent urban invasion 
flood plains: cities which have had serious flood yr. well 
cities that had unprecedented flooding only few years before, there has been 
extensive building without any immediate prospect protection. Moreover, 
certain valleys, such the Tennessee Chattanooga, the completion 
river-control works upstream has been followed further movement into the 
flood-plain where frequency has been reduced but where high flows still are 
possible. Although most Federal flood-control works are built protect against 
project flood and conceivably will one day, however infrequently, exceeded 
larger flow, there universal disposition believe that the rare flow 
will never come. This means that the number situations which catas- 
trophic disaster may follow one those rare flows new levees 
are completed. 


These generalizations are illustrated the near-by cities Boulder and 
Denver. 


FLOOD-PLAIN OCCUPANCE BOULDER 


There has been major flood Boulder since 1894, when waters from 
intense rainfall the Rocky Mountain Front area Boulder Creek inundated 
the areas the north the Creek far Water Street. (Fig. 1950, 
Congress authorized combination levee, channel improvement, and bridge 
reconstruction project protect sectors the plain that had been built 
over the years. The authorized project was accommodate flow equal 
that 1894, about 12,000 cfs. flow much 15,000 cfs would not 
impossible. The City decided not take part the authorized project 
which would have required local contribution $1,192,000 out total es- 
timated cost $1,707,000. the situation stands today with protection 
works, and with some local interests pushing for Department Agriculture 
watershed survey possible line alternative solution. 

Since 1938, the area which was flooded has been the arena for vigorous new 
construction. Single-family residential structures have increased 34%, com- 
merical structures 26%, and public structures 300%. The University Colo- 
rado has built twenty eight new multiple-unit housing structures the flood 
plain. The City has constructed its municipal building between the channel 
and the line the proposed levee. The largest bank has its new building 
astride the proposed levee line. One modern hotel under construction 
short distance below the canyon mouth, and another along the creek between 
24th and 28th streets. 

Some managers property the flood plain believe that the area never 
will flooded again. Others are taking calculated risk, arguing that the 
benefits daily use will more than offset the occasional over- 
flow, and that another flood may not come their period occupance. They 
are correct thinking that Boulder has not recently added the national bill 
flood losses. But when flow equivalent the one 1894 does come down 
the channel—perhaps next year, perhaps not our lifetime—the city will 
the front page every newspaper the country place where property 
losses exceeded several millions dollars few hours. Hopefully, the 
pounding creek waters, laden with boulders and debris, may not cause loss 
life. the experience elsewhere applies here, that unpredictable but certain 
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5.—MAP BOULDER, COLO. SHOWING APPROXIMATE FLOOD PLAIN 
AND LOCATION SOME NEW DEVELOPMENTS. 
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flood will followed organization emergency citizens committee 
obtain Federal flood control for Boulder. 


FLOOD-PLAIN OCCUPANCE DENVER 


Denver has three distinct flood situations, one which has been largely 
solved, one which daily becomes more acute, and one which course 
halting but imaginative solution. The Cherry Creek flood plain (Fig. was 
the scene damaging overflow 1933 that covered large sections resi- 
dential and commercial occupance the city. This area now protected 
Federal detention dam capable controlling storm runoff more than 
600 cfs per sq. mile. There has been argument that such protection unduly 


6.—MAP DENVER, COLO. SHOWING APPROXIMATE FLOOD PLAIN 
AND LOCATION MAJOR TRANSPORTATION LINES, 


expensive, but, regardless cost, that sector Denver seems off the 
rolls for flood-loss reporting. 

Along the main stem the South Platte River, both above and below the 
mouth Cherry Creek, there broad plain that received partial protection 
from WPA works and that has not experienced flood flows exceeding the chan- 
nel capacity 13,000 cfs recent years. The Congress, 1950, authorized 
reservoir and levee protection project control flows 85,500 cfs for this 
area. The dam was never started, and, indeed, could not started today with- 
out dislodging unit the Martin Aircraft Corporation and residential sub- 
division, both which have since invaded the proposed reservoir area. Mean- 
while, the South Platte flood plain commercial structures have increased 
26% since 1936, industrial structures 40% 50%, and residential structures 
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FLOOD-PLAIN OCCUPANCE 


slightly, notwithstanding extensive demolition for new expressway inthe val- 
ley. The damage potential unquestionably has increased rapid rate this 
lower valley area. 

little progress has been made the small tributary gulches where there 
has been persistent tendency for residential and commercial occupance 
press toward the dry stream beds, thus reducing the capacity the channels 
and exposing more property flood loss. Denver one the few cities that 
has begun program zoning, building regulation, and land acquisition 
control further encroachment such dry and start has been made 
metropolitan area planning for regulation across political 

While the Cherry Creek area flood loss has been eliminated engineer- 
ing works, the South Platte damage potential has been growing under the stimu- 
lation highway and partial flood-control construction. modest startis un- 
derway curb further encroachment the tributary gulches. 


RANGE CHOICE ADJUSTMENT FLOOD HAZARD 


The traditional approach flood losses the United States has been 
either control the flood flows engineering works bear the losses, let- 
ting the public come the aid the less prosperous flood sufferers. This 
shown the heavy public emphasis upon investment Federal flood-control 
works and the Red Cross and other emergency aid programs. 

There are, course, other possible adjustments which individual property 
managers can make. The whole range choice may listed follows: 

Bear the Losses.—A good many households and firms bear the losses when 
they occur, and few actually make financial plans so. 

Emergency Evacuation and Rescheduling.—Property may removed from 
the reach floods, and property movement and production operations may 
rescheduled avoid losses through interruption. This requires rela- 
tively accurate system flood forecasting and plan for emergency action 
when the critical forecast received. 

Prevent Flows.—In some areas practicable prevent certain flows, 
particularly the more frequent flows, land treatment and associated mea- 
sures. The experience with such measures only beginning build up. 

Elevate Land.—By land fill practicable some situations raise 
property above the level flood waters. The effect upon channel capacity 
varies according location with respect flowage and pondage areas the 
flood plain. 

Control Flows.—This the conventional engineering solution, involving 
channel improvement, levees, cut-offs, and storage detention dams. 

Change Structures.—In advance flood warning, structures may al- 
tered prevent reduce flood losses when the water rises. These 
changes may include such alterations packing machinery, bricking 
low openings, cut-off valves sewers, and rearrangement electrical cir- 
cuits. 

Change Land Use.—The use flood-plain land may changed in- 
troduce use that less susceptible flood loss. This may range from the 


“Report for Inter-county Regional Planning Commission Storm Drainage the 
Denver Metropolitan Area,” Dale Rea, 1957. 
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transfer entire town from riverine upland site, public acquisi- 
tion flood-hazard areas for recreational parking purposes. may 
guided public regulation, including zoning, building ordinances, subdivision 
regulation, and land acquisition. 

Insure.—Although insurance against flood losses not generally available 
the United States, there inactive program for Federal-state subsidized 
insurance under the Flood Insurance Act 1956, andthere are few instances 
coverage private companies where special structural adjustments have 
been made. 

Public Relief.—Through the medium the Red Cross direct grants 
and loans under Public Roads, Civil Defense, Corps Engineers and small 
business programs, public aid given individuals and local governments 
that suffer heavily from floods. 

There are concrete examples each type adjustment, and there are 
strong arguments that can made for and against each type adjustment, 
depending upon the local situation. understand how and why particular ad- 
justments are made given situation necessary look the various 


elements that enter into decisions management urban flood-plain prop- 
erties. 


ELEMENTS DECISIONS FLOOD-PLAIN USE 


with most other resource-management situations, least seven con- 
siderations seem enter into decisions flood-plain adjustments. may 
helpful examine how each one figures the urban situations which have 
been studied. 

Estimating the Resource.— There widespread ignorance the flood hazard 
and tendency minimize it. Many people building buying flood plains 
are unaware the precise hazard they are running, grossly misinterpret 
the technical estimates. This applies even places where there have been 
public plans for flood protection. man says need not worry about floods 
because 200-yr. flood occurred the year before. Federal housing insur- 
ance office says does not insure mortgages for new buildings flood plains 
but has map showing where the flood-hazard areas are located. 

Discounting Future Benefits and Costs.—It common Federal flood pre- 
vention and flood-control studies discount future benefits and costs from the 
proposed works. These cost-benefit ratios have significance Congressional 
decisions long they are below unity. There tendency group all 
programs with ratio more than unity together and pay little attention 
the means which the ratios were calculated. far private planning 
flood-plain adjustments are concerned there little evidence such discount- 
ing procedure. few managers large industrial units make cost-benefit 
calculations but most managers not. 

Harmonizing Two More the Federal plans there disposi- 
tion harmonize flood-control plans with other water uses such irrigation 
and navigation. the level local governments and private managers such 
multiple-uses tend ignored. Few attempts are made prevent harmful 
encroachments, and there little attention combining levees with highways, 
harmonizing park development with floodway improvement. One 


the hopeful moves this direction under planning supported the Urban 
Renewal Administration. 
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Projecting Future Demand.—Typically, the Federal agencies assume that 
there will little change the demand for flood-plain land while private 
overestimate the effect protection and rush the hint 
prospect some control work. 

Technological Change.— Private managers tend overestimate 
the physical benefits from any kind prevention protection work. The 
completion single-purpose power dam upstream, (and even downstream), 
the beginning watershed improvement, the clearing upstream 
channel taken assurance that floods will abated the future. Thus, 
little protection work may encourage great deal channel encroachment. 

Integrating Regional Uses.—While the Federal agencies try assess the 
consequences flood-control program upon other water uses inthe same 
adjoining basins, there attention this the part local and private 
managers, and metropolitan area plans for storm water disposal are rare. 
Even less common are efforts assess the possible use flood plains 
serving regional requirements for land and for riverine facilities. 

Setting Social Guides.—We come now one element which affects all the 
rest. Society, through its public agencies, definitely restricts some kinds 
decisions property managers and clearly encourages some other kinds 
action. Federal agencies tend guide the occupance flood plains pro- 
viding prevention and protection plans, giving public relief, issuing storm 
and flood warnings, and offering information flood occurrence. Seven 
states attempt regulate channel encroachment. few cities and counties 
exercise some regulation over flood-plain use. One Federal agency cooperates 
with state planning agencies assisting with local plans for flood-plain use.8 


BROADENING THE RANGE CHOICE 


the present circumstances the social guides the form information, 
regulation, and investment affecting flood-plain occupance tend encourage 
further encroachment upon the flood plains the same time that they lead 
heavier Federal expenditures for flood control. not entirely whimsical 
say that the national situation somewhat like local situation where investi- 
gation showed forty new houses flooded area could not economically 
protected levee works. was noted, however, that twenty additional 
houses were constructed protection then would feasible. the 
existing policies one may expect that the additional twenty houses will time 
built and that Federal largess will follow due course. 

know that the managers flood-plain properties often underestimate the 
hazard, fail discount future benefits and costs, overestimate the demand for 
their land, overestimate the effects engineering and land-treatment works, 
and ignore possible combinations flood-loss reduction with other local im- 
provements with regional plans for land and water use. also know that 
practical matter most managers who face aflood threat have 
tween only two alternatives. They can bear the losses, they can press for 
Federal project some sort. They not receive technical advice the 
possibilities emergency evacuation, land elevation, structural changes, 


“Regulating Flood-Pain Development,” Francis Murphy, Univ. Chicago 
Georgraphy Research Paper, 1958, No. 56. 
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land-use changes. They not have insurance readily available. they suffer 
unduly heavy losses they have the prospect public aid. they are energetic 
enough they may receive public protection. With few exceptions they not 
know the precise character the flood hazard, and they are not pub- 
lic curbs against any further encroachment upon the stream channel. 

are break out the present situation which floodlosses promise 
tokeep pace with flood protection for least several decades seems essen- 
tial broaden the range choice. Each manager—public private—should 
given the opportunity choose among the whole range possible adjust- 
ments, with the public agencies strictly limiting those choices which would 
cause damage others the public safety. This calls for new and fresh 
approach the problem flood loss United States. impor- 
tant beginning that direction was made the Conference called the Coun- 
cil State Governments December, 1958.9 The new approach does not re- 
quire abandonment existing public programs. requires that they 
supplemented new expanded efforts broaden the choice open 
face flood hazard. 


The chief directions which now seems necessary move are these: 


Publication and wide distribution flood hazard maps and reports for all 
important flood hazard areas, made available both property managers 
and public lending and construction agencies. 


Technical advice property managers the means and costs re- 


ducing flood losses emergency measures, structural changes, and 
land-use changes. 


cooperative Federal-state program flood insurance which pre- 
miums are proportion hazard. 


State regulation ofany further channel encroachments which would cause 
damage others the public safety, and Federal requirement that addi- 


tional Federal expenditures for flood control contingent upon such regula- 
tion. 


improved system for issuing flood warnings and for getting them into 
the hands the property manager concerned. 


Technical assistance local governments drawing plans for re- 


ducing flood losses whatever combination means may seem most feas- 
ible. 


New Attack Flood Losses,” Gilbert White, State Government, Spring, 
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APPLICATION SNOW HYDROLOGY THE COLUMBIA 


JOHNSON, ASCE! and BOYER.2—Mr. Miller has presented 
well-written, thoughtful discussion, adding several explanations for which 
the writers are grateful. They agree that all the detailed exploratory studies 
comprising major research effort should made available those inter- 
ested the general subject covered. His statement that the original papers 
are nowhere referred little misleading, that the referenced “Snow 
Hydrology” has some pages filled with references works writers 
the general field hydrology well those engaged analyzing and re- 
porting upon the results the recent laboratory measurements. Some 
these were written Mr. Miller. 

The water-balance approach forecasting volume streamflow from 

melting snow may indeed replace, time, the index methods now used. How- 
ever, much better coverage must first established such elements 
snow water equivalent, soil and snow priming, evaporation and other losses, 
and precipitation. Until such time, the index methods presently use must 
continue serve. fact, until long-term weather forecasts can made 
with improved reliability, forecasts spring snowmelt runoff will continue 
made with considerable uncertainty. 

Mr. Miller chooses reverse the order consideration four equations 
concerning the heat balance concept, going from formula containing large 
number terms one containing only two temperature measurements. For 
presentation hydrologist completely familiar with the subject, this order 
consideration might logically followed, but the writers assumed that 
some persons whose knowledge snow hydrology was limited would read the 
paper. this case, the usual and preferable method, where choice exists, 

start with the least-complicated situation and develop formulas that 
consider additional parameters, until the most complicated situation cov- 
ered. 

Mr. Miller emphasizes the sparsity radiation observations, particularly 

for long wave radiation. Several relatively inexpensive and self-contained in- 
struments for measuring and recording solar radiation have been developed, 
one English concern, based upon research the Meteorological Office. 
result, may soon possible equip stations mountainous basins 

instruments for measuring daily solar radiation; this information would 
reported the observer the using offices electric means order 
that better forecasts runoff may made. 


January, 1959, Oliver Johnson and Peter Boyer. 
Hydr. Engr., North Pacific Div., Corps Engrs., U.S. Dept.of the Army, Portland, 
Oreg. 


Hydr. Engr., Portland Dist., Corps Engrs., Dept. the Army, Portland, 
Oreg. 
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Discussion William Sammons 


WILLIAM SAMMONS, ASCE.—The writer deeply indebted the 
author for personal assistance prior publication this graphical method 
for the determination K-values.2 

April 1958, the writer was “knee-deep” the preparation 
for the May meeting the American Geophysical Union. phone call the 
author solved the major problem which dictated solution. This will re- 
ferred later Method and will presented briefly for the benefit the 
reader. Methods employed the May, 1958 paper stimulated the writing 
subsequent September, 1958. This latter paper presented theory, 
examples and K-values complied from other few new 
applications were also presented. 

The K-values presented here Table “Theoretical Log-Probability 
Critical Values,” will discussed under Method and was designed sup- 
plement available tables.4,8 

Method A.—Limited time demanded the expediency this procedure which 
the author’s suggested method (with few liberties taken). The writer be- 
lieves the K-values developed (tables the K-value were furnished the 
author and few office copies were printed) served the purpose for which they 
were intended.3 The writer does not recommend the indiscriminate use the 
unpublished table. Briefly, Method follows: 


July, 1959, Ven Chow. 
Hydraulic Engineer, U.S. Agric., Soil Conservation Service, Central Tech- 
nical Unit, Beltsville, Maryland. 

Method for Determining the Log-Normal Frequency Fac- 
tor,” Ven Chow, (Private communication with Ven Chow, April 1958. Not 
for publication). 

“Agricultural Watershed Influence Runoff Potential,” Melvin Kleen, and 
Wm. Sammons, presented the thirty-ninth annual meeting the May 5-8, 1958 
Washington, (Manuscript only). 

“Relationship Between the Original and Transformed Variates for the Logarithmic- 
Probability Law,” Wm. Sammons, (For in-Service use and distribution the SCS 
Hydrology Meeting, Jackson, Miss., October 13-17, 1958). 

the Logarithmic-Probability Law,” Kalinske, Trans. AGU, 1946, Vol. 
27, No. pp. 709-710. 

Reference Proc. Paper 2084, Ven Chow. 

Nomogram for Log-Normal Frequency Analysis,” Weiss, Trans. AGU, 
1957, Vol. 38, No. pp. 33-37. 

«Critical Values the Log-Normal Distribution,” Jack Moshman, Statis- 
tical Assoc., 1953, Vol. 48, No. 263, pp. 600-609. 

Univ. Press, 

“Fitting Generalized Log-Normal Distribution Hydrologic Data,” 
Brakensiek, Trans. AGU, 1958, Vol. 39, No. pp. 469-473. 


Re 
= 
& 
i 
4 


108 February, 1960 


The author’s procedure was followed verbatim. 

Figures were plotted similar the author’s Fig. 

Probability lines were developed for twelve “key” 
(Cy 1.1171 3.455). 


MINIMUM 


0.146 


NUMBER ITEMS EACH CLASS 


18 2) 24 27 30 O03 O89 12 15 18 24 27 30 
COEFFICIENTS VARIATION 


FIG, 1.—DISTRIBUTION THE COEFFICIENTS VARIATION, 


K-values (graphical determinations) were plotted the logarithmic 
scale versus the arithmetic scale semi-logarithmetic 
graph paper. 

K-values for nine levels probability used Table Reference 
were read from these graphs needed for use Reference 

These nine K-values for specific Cy-values were smoothed com- 
puting X/X; plotting probability paper; drawing average line 
eye; re-compufing K-values the author’s Eq. 


f-0.103 iti 
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Smoothed K-values were used calculate X-magnitudes the hydro- 
logic event. 

table approximately 100 (or Cy) values were tabulated with 
values the nine levels probability. 


Fig. presented illustrate the extent this unpublished tabula- 
tion K-values was employed. 
Method B.—K-values developed Method could have utilized IBM- 
650 computer great advantage and under different set prevailing cir- 
cumstances, would probably have been worth while official project. 


The author’s Eq. was modified follows: 
= 


where K-values for and few selected values are tabulated 
Table 


TABLE 1.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 
FOR THE STANDARD NORMAL VARIATE EXCEEDED WITH 
SPECIFIED 


Probability 


Probability 


0.000000 


Five levels probability were employed compute K-values. 
Codex Book Company, Inc., Logarithmic Normal Probability Paper No. 
31,376 was used plot X/X versus probability the 99, 84.13, 50, 
15.87 and probability levels (used smoothing). 
The five levels probability, “computed” K-values, were plotted 
and the inch, coordinate paper No. 359-10 1/2 LG, versus 
Smooth curves were drawn and computational errors located and cor- 
rected computational methods. 


-3.090232 0.10 
99.75 2.807034 0.25 2.807034 
99.6 2.652070 0.40 2.652070 
99.5 2.5758 0.50 2.5758 
2.326348 1.00 2.326348 
97.5 2.50 
+1.750868 
1.6449 1.6449 
93.5 1.514102 1.514102 
1.281552 
1.036433 1.036433 
-0.674490 +0.674490 


0.000000 


The five K-values for specific (or were smoothed the 
process described under Method step These were later plotted 


K-Value 
Probability percentage equal greater than the given variate. 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 


Coefficients Probability, percentage, equal greater than the given variate 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 


Coefficients of 


Probability, percentage, equal greater than the given variate 


Skew 


| 


3.55 
3.56 
3.57 
3.58 
3.59 


3.60 
3.65 
3.70 
3.75 
3.80 


3.85 
3.90 
3.95 
4.00 
4.05 


4.10 
4.15 
4.20 
4.25 
4.30 


4.35 
4.40 
4.45 
4.50 
4.55 


4.60 
4.65 
4.70 
4.75 


0.9220 
0.9238 
0.9256 
0.9274 
0.9292 


0.9310 
0.9399 
0.9487 
0.9574 
0.9661 


0.9747 
0.9832 
0.9916 
1.0000 


1.0165 
1.0408 
1.0488 


1.0723 
1.0800 


1.0953 


-1029 
-1104 


-1470 
-1684 


-1962 
1.2098 
1.2232 
1.2365 


1.2496 
1.2625 
1.2753 
1.2879 
1.3004 


1.3127 
1.3249 
1.3369 
1.3488 
1.3605 


1.3722 
1.3837 
1.3950 
1.4063 
1.4174 


| 1.0645 | 


.95600 

-95040 


-94215 
-93547 
.92249 


.91619 
.91004 

-90401 
-89806 


.88656 


-86477 


.85957 
-83971 


.83497 
-83026 


-82117 


.80387 


.78373 

.76875 
.76161 


.74788 


-71680 
-71106 
-69478 


-68957 
67964 


.81239 


141 
141 
140 
139 
140 


139 
686 
668 
652 
646 


630 
615 
603 
595 
581 


714 
698 


655 
639 
621 
606 


587 
574 
557 
543 
528 


521 
504 
489 
483 
469 


.72067 


-71999 
-71931 
-71800 


-71736 
-71087 


-70146 
-69544 


-68363 
-68095 
-67545 


-67270 


-66741 


-66478 


-65968 


-65219 
-64977 


.64740 


-64501 


-64035 
-63583 


-62281 
-61463 


-61067 
-60682 
-60302 
-59932 
-59567 


-58863 
-58522 
-58187 


-57537 
-57223 
-56913 


365 


355 
349 
341 
335 
326 


324 
314 
310 


-56611 302 
-56313 298 


-28766 


-28796 
-28865 


.29239 


-29291 


-29587 


-29654 
-29682 


-29811 
-29833 


-29858 
-29875 


-29934 


-29995 
-30018 


-30028 
-30026 


-29995 


-29939 
-29919 
-29902 


-66259 


.66165 


-66071 
-65995 
-65902 


-63755 
-63352 
-62580 
-62194 


-61820 
-61450 
-61065 
-60706 


-59993 
-59304 
-58956 
-58614 


-56352 
-56056 
-55156 


-52900 


-51345 
-50855 
-50371 
-49874 


-48524 
-48084 


-47228 
-46812 
-46409 
-45629 


wonr 


5.30 138 771 432 561 3.9671 
4 
5.80 129 385 503 3.9739 
5.90 128 380 490 3.9746 
6.00 126 370 484 3.9750 
6.40 120 440 3.9749 
6.70 117 430 3.9730 
7.10 111 383 3.9684 
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| Variation | 
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Probability, percentage, equal greater than the variate 


-66554 
-66102 
-65667 
-65238 
-64817 
-64407 


-64407 
-64003 
-63608 
-63226 
-62848 


-62479 
-62115 
-61757 
-61409 
-61068 


-60731 
-60401 
-60079 
-59760 
-59451 


-59143 
-58842 
-58256 


-57688 
-56875 
-56615 


-56355 
-56101 
-55851 
-55604 
-55359 


-54199 


-53978 
-53121 


-52710 
-51510 


-51733 
-51358 
-51174 


-55732 
-55453 
-55175 
-54635 


-54635 
-54373 
-54114 
-53862 
-53613 


-53368 
-53127 
-52886 


-52198 
-51975 
-51757 
-51541 
-51328 


-50913 
-50710 
-50509 
-50313 


-50116 
-49925 
-49366 


-49183 
-49006 
-48827 
-48654 
-48480 


-48313 
-47816 


-47028 
-46877 


-46578 
-46433 
-46148 


-46004 


-29833 
-29813 


-29734 
-29709 
-29679 
-29652 
-29622 


-29593 
-29562 
-29503 
-29470 


-29437 
-29405 
-29340 


-29276 
-29241 
-29208 
-29141 


-29106 
-29039 
-28968 


-28933 
-28900 
-28864 
-28830 
-28795 


-28761 
-28690 
-28622 


-28586 
-28551 


-28378 
-28310 


-45243 
-44874 
-44498 
-43787 


-43096 
-42768 
-42428 
-42115 


-41797 
-41486 
-41169 
-40569 


-40282 


-38878 
-38622 
-38361 
-38103 
-37853 


-37110 


-36180 
-35730 


-34460 


-34256 
-33863 


-33096 


-31689 


7.90 1.5026 103 395 259 328 3.9543 
8.10 1.5228 101 378 249 313 3.9499 
8.20 1.5328 100 369 245 318 3.9476 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES (Continued) 


Coefficients Probability, percentage, equal greater than the given variate 


84.13 15.87 
+ 


12.70 1.9102 -50993 181 .45321 


3.8182 


13.40 1.9593 138 334 .44429 
13.60 1.9730 137 -49462 326 .44185 
13.80 1.9866 136 -49142 320 .43946 
14.00 2.0000 134 -48831 .43712 
14.20 2.0133 133 -48526 305 .43482 
14.40 2.0264 131 -48228 298 .43258 


14.40 2.0264 131 -48228 298 .43258 
14.60 2.0394 130 -47938 290 .43037 
14.80 2.0522 128 -47654 284 .42822 
15.00 2.0650 128 -47373 281 .42609 
15.20 2.077 126 -47100 273 .42401 


15.40 2.0900 124 -46835 265 .42198 
15.60 2.1020 120 -46580 255 .42003 
15.80 2.1146 126 -46316 264 .41801 
16.00 2.1268 122 -46063 253 .41606 
16.20 2.1388 120 -45816 247 .41416 


16.40 2.1507 119 -45575 241 .41230 
16.60 2.1625 118 -45337 238 .41046 
16.80 2.1742 117 233 .40865 
17.00 2.1858 116 -44876 228 .40688 
17.20 2.1973 115 -44652 224 .40513 


17.40 2.2086 113 -44433 219 .40343 
17.60 2.2199 113 -44216 217 .40173 
17.80 2.2311 112 -44004 212 .40007 
18.00 2.2422 111 -43795 209 .39844 


18.40 2.2642 109 -43387 201 .39523 
18.60 2.2750 108 -43190 197 .39367 
18.80 2.2858 108 -42994 196 .39213 
19.00 2.2965 107 -42801 193 .39061 
19.20 2.3070 105 -42614 187 .38913 


19.40 2.3175 105 -42429 185 .38766 
19.60 2.3285 110 -42236 193 .38613 
20.00 2.3486 103 177 .38336 
21.00 2.3988 502 -41043 845 .37661 


22.00 2.4473 485 785 .37029 
24.00 2.5397 455 -38841 686 .35878 


247 .27836 .30411 
244 .27771 .30103 
239 .27705 .29804 
234 .27639 .29513 
230 .27574 .29231 
224 .27509 .28956 


224 .27509 .28956 
221 .27447 .28678 
215 .27382 .28422 
213 .27321 .28154 


-25600 214 3.6815 


675 .25647 251 .22205 795 3.5777 263 


593 .25173 233 .20781 684 3.5286 250 


306 3.7967 
308 3.7905 
299 3.7844 
291 3.7783 
282 3.7721 
275 3.7659 


275 3.7659 
278 3.7598 
256 3.7536 
268 3.7475 


wo 
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probability paper, noted under step above, check linearity. 
Forty (or Cy) values proportionally spaced throughout Table con- 
stituted the sample tested. 

First differences, are shown Table aid prospective 
users. 


Needless say, the five computed K-values for each entry will al- 
low further probability levels obtained the author’s Fig. method 
utilizing tabulated K-values for specified values and solving the author’s 
Eq. (If publication space were available the columns, lines and 
pages tabulations the individual steps could made available). Table 
does not contain additional forty-eight computed which have 
been computed date for specific problems and satisfy the curiosity the 
writer. brief summary similar work the writer noted above 
available.9 

Precaution Necessary.—As pointed out detail and K-val- 
ues may used with the original variate with log-transformation vari- 
ate. The latter more efficient, statistically speaking, than the scalar trans- 
formation the original variate. Alden Foster method for 
making theoretical corrections Similar work the writer also 
supports this generalized approach. 


Discussion “The Log-Probability Law and Its Engineering Applications,” 
Alden Foster, 1959, Proc. No. 709, Vol. 81, pp. 14-19. 
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ROLL WAVES AND SLUG FLOWS INCLINED OPEN 


Discussion Hugh Taylor and John Kennedy 


HUGH TAYLOR,! A.M. ASCE and JOHN A.M. ASCE.— 
The author has explored obscure byway hydraulics and has presented 
interesting body experimental data resulting from what was obviously 
painstaking and thorough laboratory investigation. While knowledge valua- 
ble per se, the fact that the phenomena dealt with here are little known 
confers additional value upon the results, they offer the engineer means 
broadening his intuition area not likely have experienced di- 
rectly. this respect, then, the author has done the profession definite 
service. 

regretted, however, that the author’s treatment his results 
leads one conclude that they are more directly useful than they actually 
are. This misunderstanding seems have two different sources. First, the 
choice nomenclature unfortunate that the term waves” already 
has generally accepted meaning, not precisely defined one. spite 
this, however, the author has chosen apply the term exclusively insta- 
bilities laminar flow, phenomenon which predominately academic 
interest the hydraulic engineer. Furthermore, the term defined the 
paper exclude the meanings already well established the literature, 
notably the writings and Dressler.5 Second, the confusion 
enhanced the reference impact and freeboard problems channel de- 
sign caused instabilities flow. Certainly, the laminar instabilities the 
author calls waves” are not this magnitude, they not occur 
Reynolds numbers greater than 415. 

The author has likewise limited himself choosing include within his 
term “slug flow” only those unsteady flows which have arisen more less 
directly from laminar flow. This may well reservation without reason, 
the properties the limiting wave are not influenced the presence 
absence prior laminar flow. Unless this true, however, the similarity 
these waves the roll waves reported those discussed 
Dressler® questionable. seems improbable that the flow with roll waves 
which Dressler attributes Cornish, the flow with roll waves observed 
Vanoni and the writers the Santa Anita Flood Control Channel, near 
Arcadia, Calif., could ever have been laminar. This not object the 
scope the paper, but merely question the relation between the actual 
scope and that which suggested the title and introduction. 


July, 1959, Paul Mayer. 
Dept. Civ. Engrg., Calif. Inst. Tech., Pasadena, Calif. 
Dept. Civ. Engrg., Calif. Inst. Tech., Pasadena, Calif. 
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Purpose.—The purpose this discussion three-fold and may sug- 
gested passage from the paper itself: “Rigorous mathematical treatment 
complex hydraulic phenomena encounters great difficulties. Convenient ap- 
proximations often lead equations questionable veracity. The background 
experimental reality provides important insight into physical phenomena.” 
Specifically, the writers’ intent call attention some theoretical argu- 
ments recent English paper, enter objection what the author pre- 
sents the theory roll waves” and the presentation and analysis his 
experimental results, and mention some recently reported experiments 
which supplement the findings the author. The objections will entered 
first. 

The first objection that the theoretical development weak. oneassumes 
irrotational flow and surface disturbances small amplitude with respect 
their wave length, then the linear analysis Airy2 follows. well known 
that this analysis the celerity, wave given 


where “kinematic surface tension,” called analogy with kine- 


matic viscosity, anda wave number 


Under the usual simplifications, tanh taken equal whena 
the former case one speaks shallow water waves, the expression for the 
celerity which has useful minimum; that is, 


from which 


The latter vanishes only for which requires infinite wave length 
disturbance all. the two simplifications, only the latter one, 
where one speaks deep water waves, that one finds minimum clerity and 
this case certainly irrelevant the problem hand. The analysis Eqs. 
(1) through (6) for the deep water case, and since the author makes almost 
use his discussion wave celerity, perhaps would have been better 
omit completely. 

The author’s criterion for the formation progressive waves very hard 
accept the present situation. the progressive waves associated with 
“roll waves” and “slug flow” are treated discontinuous surge fronts, and 
the velocity assumed constant with depth, the celerity the waves can 
shown 


where the velocity the wave relative the undisturbed flow, repre- 


sents the depth undisturbed flow, and the depth upstream from the 
disturbance. the disturbance vanishingly small, that is, then 
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the celerity becomes that given the author Eq. 11, The same 
limiting value given Eq. for shallow water waves surface tension 
neglected and the disturbance very small compared the wave length 
and the depth, that is, the deep water case. 

For Eq. valid, the velocity must not vary with depth, upstream and 
downstream from the disturbance. some other velocity profile, such the 
parabolic one used the author, exists upstream and downstream from the 


disturbance, the wave celerity further modified, was shown 


where the amplitude the disturbance, represents the Froude Num- 


and denotes the velocity the disturbance relative the average velocity 
the undisturbed flow. 

lar velocity profile). Acomparison Figs. and indicates that the ex- 
periments reported the author, the depth flow and the wave heights were 
the same order magnitude. Noticing also that the Froude numbers were 
quite high, and recalling that the momentum coefficient has the value 1.20 for 
parabolic velocity distribution, one must conclude that the use 
the wave velocity relative the average velocity flow incorrect. 

Actually, taking the Froude number 1.7, the amplitude the disturb- 
ance half the undisturbed depth, parabolic velocity distribution, the 
resulting disturbance velocity only about half that given the limiting 
case. Similarly, the celerity oscillatory waves (Eq. greatly modified 
the velocity profile departs appreciably from rectangular. Thus Eq. and 
the criteria stated Eqs. 12a and must rejected. 

The kinematic argument Lighthill and properly predicts 
celerity very long, “kinematic” waves (in which dynamic considerations 
are unimportant) equal twice the mean stream velocity. But convert this 
the statement that the Froude number equals equivalent stating that 
the wave celerity given which was shown above incorrect for 
the problem hand. and Whitham themselves said: “It must re- 
marked that the result has been deduced under the assumption that the 
disturbance small. More generally,..., the value which must 
exceeded bore constant strength maintained, depends the 
strength; this value always less than and tends the strength 
approaches zero.” 

The second objection entered that the author’s analysis the ex- 
perimental data somewhat misleading. Where dimensionless parameters 
are used, generally desirable use the standard ones for which the 
reader already has intuitive “feel,” resist the strong temptation 
invent new ones unless the gain greater simplicity the new ones out- 
weighs familiarity with the established ones. The parameter 


presented Eq. which calls the critical flow number actually 
modified Weber number: 


sur Théorie des Eaux Courantes,” Mémoires Divers Savants 
mie des Sciences, 1877, Vol. 23. 
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This results from the tacit assumption that the Froude number unity the 
development the critical flow number Eq. 29. 


The critical depth number, Np, which given Eq. 31, can also con- 
sidered shear velocity Reynolds number, 


Np = (5 ) = (v) 


where shear velocity DS. The Darcy-Weisbach friction factor, 
defined 


and can expressed 

from which 


the laminar regime, where parabolic velocity profile applicable, Eqs. 
10a and (vi) yield 


the present experiments, single fluid was used and thus the only vari- 
ables the critical flow number are Thus Fig. 16, which 
plotted against essentially relates and and would expected 
give linear relation the laminar region logarithmic plot. With 
suitable manipulation variables, Fig. can put the form the usual 
pipe friction diagram (except for the constants involving and which 
found any elementary text fluid For flow this type, the 
laminar region ends approximately 600 1,000 (although laminar 
flow can exist Reynolds numbers 2,500 under appropriate conditions) 
which corresponds 42.5 54.8. Fig. 16, this about the region 
which the straight line relation ends, would expected. also the up- 
per limit the Reynolds number for which roll waves were observed and this 


confirms the conclusion that roll waves” are confined the laminar regime 
flow. 


«Fluid Mechanics,” Daugherty, and Ingersoll, McGraw-Hill, 1955, 
182. 


g1/6 
Substituting this value Eq. viii gives 
t 
a 
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similar manner, Figs. 12, 13, 14, 15, and are also essentially pipe- 


friction diagrams for this type flow. This becomes more apparent when 
considered that 


Indeed, use made Eq. the result for the laminar regime 


Ne = VS/3 Rp 


which exactly the relation the author presents Eq. and Fig. 18. The 
runs conducted outside the laminar regime are the transition regime 
from laminar turbulent flow, thus the friction factor cannot obtained di- 
rectly from pipe friction diagram, since the relation between and 
this region not well defined. Actually, depends the manner which the 
transition occurs, which turn depends the manner which the experi- 
ments were conducted. Had the experiments progressed into the fully turbu- 
lent regime, the relation between and could obtained directly from 
pipe friction diagram. Thus appears that this phase the author’s analy- 
sis consists primarily presenting the well known pipe friction diagram 
several alternate forms. 

The author’s statement that the minimum slope for roll wave” formation 
was not actually verified his experiments since runs were made 
with slopes the neighborhood Actually Figs. 16, 17, and would 
indicate that better criterion for the formation roll waves” and slug flow 
defined the author would that roll waves” are associated with lami- 
nar flow, while slug flow occurs when the flow the transitional regime 
between laminar and turbulent flow. The writers also question the applica- 
bility the statement: “All Froude numbers leading the formation slug 
flows were above the theoretical minimum for slug 
flow formation mentioned above would indicate that slug flow could result even 
the flow were transitional turbulent. experiment was con- 
ducted this range either prove disprove the necessity turbulence 
for slug flow formation regardless the size the Froude number. Such 
experiment would very enlightening. 

Comparison with Other Recent Papers.—There are two papers (by 
Binnie20 and which appeared about the same time* 


“Instability Slightly Inclined Water Channel,” Binnie, Fluid Mech., 
1959, vol. pp. 561-570. 

“Wave Formation Laminar Flow Down Inclined Plane,” Fluid Mech., 

Benjamin, 1957, vol. pp. 554-574; Numerical corrections ibid. 1958, vol. 657. 

The precedence somewhat confusing. The theoretical paper Benjamin 

was published before the present paper (but after the thesis upon which based), 

Similarly, the experimental results Binnie appeared after the oral presentation the 

present paper but before was published. aside, might suggested that the 

date submission paper shown case questions priority ever should arise. 

This common procedure many technical journals. 
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the present one, and which are relevant this discussion. They are men- 
tioned some length here because they appeared relatively new publica- 
tion, the Journal Fluid Mechanics, which published England, and may 
not familiar some readers. The two papers are essentially contempo- 
rary with the author’s; interesting see they are agreement. 

The English papers supplement the present one that (a) they are con- 
cerned with the lower threshold waves laminar flow, (b) Binnie reports 
having explored the spectrum Reynolds numbers from before occurrence 
waves into the turbulent patch “slug flow” condition, and (c) Benjamin 
has developed simplified stability analysis, not dependent upon the Airy as- 
sumptions, which Binnie has applied laminar flow slightly inclined open 
channels. 

brief review Benjamin’s stability analysis presented Binnie, 
adapted for small slopes. essentially follows: Define the following di- 
mensionless quantities: 


Celerity Cj, and both being real 

Time D/V. 


Now one imposes small disturbance exp i(X T)} onthe 
laminar flow represented Eqs. through 10, there results equation 
known the Orr-Somerfeld This equationcan solved power 
series which converges fairly rapidly for small although the alge- 
bra quite tedious. For very small waves one finds (thus confirming 
Lighthill and Whitham under the given assumptions), and 


where the inclination the channel. Within the small-wave assumption, 
all this equation can show the relation between Rp, and under which 
the waves would start grow exponentially. The condition for stability that 
that the disturbance will have part which grows exponentially 
with time. Thus, for instability, must have 


interest the wave number, for the fastest growing instability. This 


this last result, together with the primary flow Eqs. and substituted 
into Eq. xiii, the maximum value seen 


«The Theory Hydrodynamic Stability,” Lin, Cambridge Univ. Press, 


ag 
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seen that the minimum for which instability can exist obtained 
setting max =0. This gives 


Because the linearizations involved its derivation, the development 
mentioned above cannot lead conclusions concerning fully-developed 


TABLE 1.—THRESHOLD INSTABILITY. COMPARISON CALCULATED 
WITH OBSERVED RESULTS, TERMS DEPTH REYNOLDS NUMBERS 


(degrees) 


Note: given cot Eq. xv. Rp" the threshold reported Binnie. 


“roll waves.” However, can indicate method attack, and can also 
used give lower bound the values below which unstable waves, 
and the roll waves” arising from occur. For this purpose the 
Reynolds numbers corresponding neutral stability (that is, were 
calculated, and the values superimposed the author’s Fig. ob- 
served thresholds are also shown, for 2°, and (See Fig. 27.) 
These values are well below the author’s observed values and therefore can 
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Depth, 


0.05 


FIG, 27.—COMPARISON 
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Mean Velocity, fps 


Lower threshold 
Calculated 
Observed 


MAYER’S RESULTS WITH LOWER THRESHOLDS 
REPORTED BINNIE FIG, 17). 
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FLOOD CONTROL ASPECTS CAUCA VALLEY 


Discussion Gordon Williams 


GORDON WILLIAMS, planning and final design water 
resources developments should preceded thorough study the basic 
meteorologic and climatologic characteristics the general region. Sucha 
study should first lead classification the climate according type 
along the lines suggested Képpen, Thornwaite and others. evident that 
such study was made for the Cauca Valley, and although not specifically 
mentioned, falls into “tropical wet and dry” savanna” 
classification. Such climates extend along the entire western coast Central 
America and into Mexico. 

The purpose classifying the climate type determine gen- 
eral way the potentialities the climate, including the possible occurrence 
different storm types, the probable severity and duration droughts, and the 
possibility killing frosts and damaging winds. For example, the tropical 
wet and dry climate can expected have total annual precipitation which 
will enough supply the annual evapo-transpiration expected from 
consideration the regional temperature and humidity, but the distribution 
the precipitation month month will vary widely and there will seasonal 
deficiencies which will result crop failures without irrigation. The Cauca 
Valley definitely has this type climate and, therefore, irrigation works are 
necessary develop the agricultural potential the region. 

The paper reveals that the storm rainfall potentials the region are low, 
relative many other parts the world, because the absence hurri- 
canes cylconic disturbances associated with cold fronts. measure 
this characteristic, pointed out that the maximum 24-hr rainfall has 
never exceeded the average annual rainfall. This contrast the 
east coast the United States where 24-hr storm rainfalls equal 
20% the average annual rainfall have occurred over areas several thou- 
sand square miles. Presumably most the precipitation the Cauca Valley 
results from convective cells, which are characteristics the doldrum re- 
gion. The peak this convective activity occurs Colombia, the northern 
hemisphere’s winter, spring late fall seasons, when the solar energy the 
vicinity the equator greatest. The seasonal rainfall patterns result from 
alternate north and south migration the centers convective action 
unison with the causative solar influences. mention made monsoon 
conditions, but undoubtedly the mountain barriers the west would prevent 


the development any strong moisture-bearing monsoon winds from the 
Pacific Ocean. 


September, 1959, Phillip Kirpich and Carlos Ospina. 
Prof. Hydr. Engrg., Mass. Inst. Tech., Cambridge, Mass. 
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The method determining the spillway design storm and flood unusual 
and dictated bya lack data and the meteorologic characteristics the 
storm rainfall. The method has considerable justification but some the de- 
tails the analysis might questioned. The attempt evaluate the 1,000- 
frequency for the 5-day precipitation from the 18-yr record mathe- 
matically unrealistic. The rainfall record short that the slope the 
frequency curve very uncertain and values taken from long extrapolation 
may not necessarily conservative. spillway design storms the object 
conservative and resort the factor safety principle may the 
only rational approach. The method expressing the rainfall depths 
percentage the annual rainfall could adopted for all durations from 
days. Then comparisons could made with corresponding percentages that 
have been experienced similar climates other parts the world, where 
the rainfall records are much longer. the climate river basin has been 
correctly classified, the writer believes that statistical analogies may 
made even from one continent another. 

The procedure plotting the known 5-day rainfall totals against the max- 
imum daily flow practical and expedient view the lack records and 
the non uniform character the rainfall. However, the position the en- 
velope curve subject individual judgment. questionable whether 
straight line semi-log paper proper. the plot transferred 
arithmetic scale, the envelope curve will found concave upward the 
general form parabola with vertical axis. the infiltration losses for 
the storms are essentially constant, the envelope curve certainly violates the 
unity-hydrograph theory which shows that hydrograph peaks vary about the 
first power the volume rainfall excess. natural storage, either 
channel, valley lake type, present there will tendency for the en- 
velope curve concave downward arithmetic scale. the basis 
the above discussion, and assuming that 186 represents the correct fac- 
tor safety for the rainfall difficult justify daily flow excess 
2,500 per sec. The computed peak discharge 6,000 per sec with 
Myers rating 4,550, appears place the river climatic province 
which obviously not. one schooled the United States method de- 
termining the maximum possible flood for spillways for major dams, the 
writer would the last one deplore the use factors safety spillway 
design, but, this case, appears that the planners the project not 
have faith their own very exhaustive meteorologic and hydrologic appraisal 
the region. 

The adopted plan for flood protection and drainage the Aguablanca Proj- 
ect, benefits from the application hydrologic and hydraulic ingeniuity 
well from the sound economic analysis. The intercepting canal the west, 
which carries the three mountain streams around the protected area, good 
solution which undoubtedly aided favorable land slopes southerly di- 
rection. The use excavated pond reduce pumping capacities good 
practice which has been applied the drainage leveed areas the United 
States, and, some cases, reduce the size storm sewers ordinary 
urban drainage problems. 
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FRICTION FACTORS CORRUGATED METAL 


pipe flow accurate experimental dataon loss head corrugated 
pipes large size. regretted that the experiments could not ex- 
tended cover larger range Reynolds numbers, show the trend 
the vs. curves before and after the maximum point observed, and 
particular show the friction coefficient becomes constant higher 
Reynolds numbers. 

With reference Fig. the original paper, interesting note 
that, the recommended values friction coefficient (corresponding the 
maximum observed values) are considered representing the values for 
rough turbulent flow, the equivalent pipe roughness, based these values, 
nearly the same for the three pipes and agrees with that computed for the 
smaller pipes reported. Taking Nikuradse’s equation: 


for the definition the equivalent roughness the values computed for the 
various pipes reported are given Table 


TABLE 


Nominal 
diameter, 
feet 


September, 1959, Marvin Webster and Laurence Metcalf. 
Hydr. Lab., Technion, Haifa, Israel. 


0.054 4.303 38.38 0,18 
3.00 0.077 3.604 17,09 0.18 
1.5 0.102 3,131 9.92 0,15 
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HEAD WEIGHTED VALUE 
Discussion Manuel Benson 


MANUEL ASCE.—The revised computational procedure 
proposed the authors apparently motivated the difficulty handling 
the extremely large numbers encountered computing the velocity-head co- 
efficient This difficulty easily resolved coding process. 

Table shows the slide-rule computation for cross-section which 
has been subdivided into five parts the basis hydraulic radius and rough- 
ness. Cols. and show the areas (a) and conveyances (k) the separate 
subsections. 

Lines are drawn arbitrarily through the and columns, effect shifting 
the decimal point, that k/a-values are, general, small numbers greater 
than 

Col. shows the slide-rule computations All the significant fig- 
ures Cols. and are used for this computation. Cols. and are 
rapid mental computations based only the figures the left the lines 
the and columns. These computations are made order locate the 
decimal points the slide-rule computations The correctly po- 
sitioned figures are shown Col. Cols. through need not recorded; 
they are shown only illustrate the mental processes involved locating the 
decimal points. 

The values Col. are coded values and not represent the actual val- 


ues the cubed conveyance divided the squared area. However, com- 
puted 


(k3/a2) 
K3/A2 
dimensionless and correct regardless the coding used. 

the geometric mean the conveyances both ends reach used 


1 
the discharge formula Ky, possible develop formula for 
which makes unnecessary trial-and-error procedure such the authors 
use. The formula is: 


Hydr. Engr., Research, Geol. Survey, Washington 25, 
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which the discharge cubic feet per second; and are the total 
conveyances section1 (upstream) and (downstream) respectively— 


the length reach, feet, between sections and denotes the 
water-surface drop, feet, between sections and the constant 


gravity, 32.2 fps2; and represent the velocity-head coefficients 
sections and and the coefficient representing the percentage 


TABLE 
Approximate values of----- 
(1) (2) (3) (4) (5) (6) (7) (8) 


recovery velocity head. For expanding reaches, commonly as- 
sumed 0.5; for contracting reaches, assumed 1.0 0.9. 
reach contracting the bracketed quantity the denominator positive, 

This method, like that the author’s, permits computation velocities 
and discharges each subdivision. necessary only proportion the 
total discharge the basis the relative conveyance each subdivision. 

Using the base data the author’s sample computations (areas, convey- 
ances, and length reach), the writer took min compute the two values 
and discharge 10,700 cfs means the suggested formula for di- 
rect computation discharge. 

The U.S. Geological Survey has developed similar formulas for multiple- 
section reaches covering rather complex conditions, such alternate con- 
traction and expansion abrupt contractions bridge opening. These 
formulas provide direct solutions and eliminate trial-and-error procedures. 
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REVISED COMPUTATION VELOCITY HEAD WEIGHTED VALUE* 


Discussion Marcel Bitoun 


MARCEL BITOUN,! ASCE.—Just keep the record straight, Ber- 
nouilli’s theorem, the Appendix Messrs. Lara Schroeder’s paper, 
not correctly applied. The eddy loss always subtractive term: 


where the signs indicate that only the absolute value the difference 
hy, must considered, irrespective its algebraic sign. 


Then, for contracting reach, 0.10: 


Chf. Design Branch, Div. Flood Control, Pa. Dept. Forests Waters, Harris- 
burg, Pa. 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL 


Discussion Blench 


ample the extreme swing away from formulas based mainly theoretical 
speculations and idealizations, empirical fitting curves that, apart from 
avoiding dimensional absurdity, seem unrelated any dynamical con- 
siderations; however, they seem biased favor the belief that 
Manning’s formula the proper type, for most alluvial channels with 
dunes. neat fit has been obtained miscellany slightly incongruous 
data the patient selection sufficient indices non-dimensional groups 
variables; the various graphs, combined with the references sources 
data, should most useful laboratory research workers. 

However, for various reasons, the writer would not use the various curves 
practice. First, the data are mainly from laboratory flumes which b/v 
lies between about 105 and 106 and the breadth-to-depth ratio, b/d, prob- 
ably round 3.0 less; the field b/v usually lies between 106 and 108 
with b/d from about 30. Some the data (Gilbert, Figs. 3,4) are for un- 
naturally uniformized bed material; some are for bed-conditions that had not 
become steady (Ref. 29); some are for mixed sediment with indeterminate but 
probably significant amounts suspension (Refs. 23,29). The writer has al- 
ready pointed out2 that extrapolation such data leads results differing 
greatly from field experience certain respects. Second, the curves and in- 
dices seem have been adjusted favor Manning’s formula there 
large body field evidence that Manning’s index 2/3 should re- 
placed 3/4 overcome the difficulty that his function hydraulic 
radius; the writer has experience canal system where, for apparently 
identical bed, had specified values from 0.018 for the main channels 
(10,000 cfs) 0.0225 for the small laterals (10 cfs)—these values, based 
experience, correspond very closely the indicial adjustment just men- 
tioned. Third, cognizance taken bed-load charge which important 
the Gilbert experiments, many river problems, and occasional special 
canal ones; for example, one test showed that altered from about 
0.021 down 0.017 for given channel breadth and 0.25 sand con- 
sequence changing bed-load charge from about zero 0.04% weight, 
and increasing the slope permit carriageof the load. (Note that Manning’s 
may keep fairly constant for various discharges one site one channel 


November, 1959, Liu and Hwang. 
Prof. Civ. Engrg., Univ. Alberta, and Cons. Engr. 


Design Alluvial Channels,” Ning Chien, Proc, ASCE Separate 611, 
Discussion Blench. 
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because the effects varying charge and side shear varying stages coun- 
teract each other some extent and the incorrect index, 2/3, acts reduce 
the variation these are not the conditions the problem just mentioned.) 
Fourth, there are simple and dynamically satisfactory formulas for mobile- 
bed channels,4 based enormous collection field data, that solve all 
the practical problems the present paper, plus others involving bed-load 
charge. 

Although the so-called regime have the practical advantages 
simplicity and good accuracy any the field—since they were dis- 
covered analyzing major field data—they also have the fundamental ad- 
vantage being consistent with the friction-factor diagram rigid-boundary 
hydraulics, which they generalize. fact, the regime flow formula can 


written for study (it written differently and very simply for design pur- 
poses) as: 


which generalized form Blasius’ equation, should since the 
boundary formed from the water-sediment complex just the so-called 
“smooth” rigid pipe boundary formed from the pure water. bed-load 
charge and constant.) Moreover, when consider that bed-factor 
can defined V2/d and side-factor V3/b, the above equa- 
tion can converted into: 


which measures equivalent roughness height and equal 


1 

should since the bed covered fact, regime channels are 
rough and smooth the same time, drawing attention the artificiality 
these adjectives, used hydraulics, and the fact that there should 
universal flow formula which breaks into special forms for special boundary 
phases merely expressing the equivalent roughness height properly; this 
consistent with the belief universal velocity distribution, independent 
boundary phase, held authorities rigid-boundary hydraulics. Revert- 
ing rigid-boundary speculations, one them converts Eq. into the form 
the Blasius equation replacing laminar film uni- 
versal flow formula exists, and Eq. correct, then Manning for rigid 
boundary conditions must amended have 3/4 index instead 2/3. 

Statistically the writer would like see confidence bands all charts 
draw attention the amount variation that might reasonably expected. 
Although this necessary statistical refinement essential for proper under- 
standing—its omission comparable with describing the rainfall place 
its long-term average without any indication the scatter about the aver- 
age—it seldom used; few hydrologists have begun expound its value 
order prevent the usual misinterpretation the significance flood fre- 
quency curves. theabsence confidence bands the writer can only give his 
opinion, based experience plotting hydraulic data, that quite different 
looking set indices would produce just good fitting curves and, par- 


ticular, would set based the belief that Manning’s 2/3 ought re- 
placed 3/4. 
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DIVISION ACTIVITIES 
HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 
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THE PURPOSE THE HYDRAULICS DIVISION 


The purpose the Hydraulics Division has been revised include its 
newly assigned responsibilities the field water resources. The revised 
statement purpose follows: 


“The advancement and dissemination scientific and engineering 
knowledge all branches hydraulics, hydrology, hydraulic engineer- 
ing and water resources. particular this shall embrace meteorology 
and hydrology asthe sciences dealing with the occurrence the 
atmosphere, the earth surface and the ground, fluid mechanics for 
the understanding all flow phenomena, applied hydraulics for the de- 
sign and planning hydraulic structures and comprehensive sys- 
tems, and those social, and administrative aspects basic the 
conservation and utilization water essential natural resource.” 


The structure the Hydraulics Division has been revised include these 
new responsibilities. The Committee Hydrology was abolished and three 
other committees have been established. The committees established and 
their membership are follows: 


Committee Hydrometeorology 


Purpose: make studies precipitation phenomena, including their distri- 
bution space and time, cyclical variations ranging from floods droughts, 
both natural occurrences and artificial controls modification; forecasting; 
including data acquisition, compilation, analysis, and distribution. 


Linsley, Chairman 
Hiatt 

Arno Lenz 

Waananen 


Committee Surface-Water Hydrology 


Purpose: study all aspects precipitation and runoff from the time rain- 
water other precipitation touches the ground; study the natural occur- 


Note.—No. 1960-5 part the copyrighted Journal the Hydraulics Division, Pro- 
ceedings the American Society Civil Engineers, Vol. 86, February, 1960. 
Copyright 1960 the American Society Civil Engineers. 
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rences floods and droughts, snowmelt runoff, the establishment design 
floods, water losses through evaporation, transpiration, and infiltration, de- 
tention water snow deposits, etc. 


Alden Foster, Chairman 
Ackermann 

Koelzer 

McDonald 

Jens 


Task Force Spillway Design Floods 


Harvey Banks, Chairman 
Ackermann 
Paul Baumann 
Walter Brown 
Elliott Childs 
Alden Foster 
Raymond Hill 
Julian Hinds 

John Keysor 
Harold Ogrosky 
Riter 
Franklin Snyder 
David Todd 
Gordon William 


Committee Groundwater Hydrology 


Purpose: study the natural occurrence groundwater, the effect con- 
servation quantity the detention water the soil, improvement the 
groundwater balance through natural processes, artificial supplement and re- 
charging, the return groundwater stream flow, etc. 


Todd, Chairman 
Paul Baumann 

Harris McDonald 

Guyton 


The fourth committee, “Water Resources Planning,” listed page the 
December 1959 Newsletter has not yet been approved. 

The changes the Division structure and its responsibilities are direct 
result the Society’s Water Resources Coordinating Committee, active since 
1958. Mr. Tiffany, member this committee, ably represented the 
Hydraulics Division. Herbert Riesbol, Denver, Colorado, has been appointed 
Water Resources Coordinator coordinate the work the new committees. 


New Task Forces Organized Other Committees 


Committee Hydraulic Structures 
Task Force the Hydraulics Culverts 


Charles Jones, Chairman 
Ralph Weaver 

Chow 

Morris 
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Task Force the Hydraulics Bridges 


Tracy, Chairman 
John Herbich 

Delleur 

Liu 


The Task Force Water Measurement was dissolved. 


Committee Flood Control 


The Task Force Flood Plain Zoning has been renamed The Task Force 
Flood Plain Regulations. The members this task force are: 


James Goddard, Chairman 
George Shanklin 

Frank Kochis 

William Wise 


Committee Sedimentation 


The task force Erosion Cohesive Materials has been organized with 
the following membership: 


Irving Dunn, Chairman 
Enger 
Ernest Smerdon 
Frank Masch 


Member Honored for Research Work 


Dr. Norman Brooks, Associate Professor the California 
Technology, Pasadena, California, was recently notified that the Board Di- 
rection had approved the recommendation its Committee Research for 
one the Research Prizes the Society recognition his accomplish- 
ments. 

Professor Brooks was recommended the Research Committee the 
Hydraulics Division for Research Prize the Society Committee Re- 
search. The Prize consists appropriate certificate award and 
check for $100. 

The five Research Prizes the Society will formally awarded the 
Research Luncheon held the ASCE Convention, New Orleans, March 
7-11, 1960. 


Executive Committee Meets 


The Executive Committee met Friday and Saturday, November and 
1959, the Palmer House, Chicago, Illinois. Those present were: Dr. Arthur 
Chairman; Mr. Maurive Dickinson, Vice Chairman; Professor Carl 
Kindsvater; Mr. Eugene Fortson, Jr.; Mr. Harold Martin, Secretary; 
and Mr. Samuel Baxter, Board Direction contact member. 

The following are quoted from the minutes the Executive Meeting: 


Annual Convention, Washington, Mr. Fortson reported that all the 
sessions the Hydraulics Division were well attended; some started slowly, 
but session had less than attendance. 
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Publications. Mr. Baxter gave brief report action taken recently 
the Board Direction. Some the highlights are: The Society staff has been 
instructed improve appearance journals; the variousdivisions will 
quired raise standards insure that only deserving papers are printed; 
standards are understood the level present Transactions and the 
Society will publish annually bound volume proceedings papers. (See 
provement Publications” this issue the Newsletter.) 

Research. The minutes the third meeting the Hydraulics Research 
Committee were reviewed briefly; the committee’s work was commended for 
its influence attitude toward research the Hydraulics Division. 

ASCE Committee Society Fellowships. was noted that Mr.E. Fort- 


son, Jr., the incoming member, replacing Vennard; Straub the 
new Chairman. 


Improvement Publications 


the Annual Convention Washington, C., October 19-20, 1959, the 
Board Direction authorized its Committee Publications improve the 
content and appearance the Division Journals by: 


Numbering pages consecutively, alining the side edges, re-copying 
mathematics and tables well text, and re-touching author’s illustration 
copy much possible—for uniform and neat appearance; 

Accepting only the papers that can recommended for inclusion the 
annual Transactions, tentatively comparable Transactions grade above 
70; and 

Assembling, the annual bound volume Papers” (called 
“Transactions”), all Journal papers that have attracted constructive discussion. 


The following guiding memoranda would helpful Division officers, re- 


viewers, and authors. Copies may obtained from the New York office upon 
request. 


“Publications Procedures” 
“Standards for Journal Papers” 
“Discussion Rules” 


NEW ORLEANS CONVENTION 
MARCH 7-11, 1960 


Session Hydrology: 


Presiding: Eugene Fortson, Jr., Member, Executive Committee, and Pro- 
fessor David Todd, Member, Former Hydrology Committee, 


and Chairman, Ground Water Hydrology Committee, Hydraulics 
Division. 


Appraisal the Stream-Gauging Network the United States 


John McCall, District Engineer, Geological Survey, Tren- 
ton, New Jersey. 


Seepage Through Layered Anisotropic Porous Media 
David Todd, Associate Professor Civil Engineering, and 
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Jacob Bear, Graduate Research Engineer, University California, 
Berkeley. 


Pressure Distribution Partially Saturated Porous Media During 
Steady Flow Immiscible Fluids 


Verne Scott, Associate Professor Irrigation, University 
California, Davis. 


New Emphasis Water Resources ASCE 
Carl Kindsvater, Past Chairman, Executive Committee, Hydrau- 
lics Division. 

Session Hydrodynamics: 


Presiding: Eugene Fortson, Jr., Member, Executive Committee, and Fred 


Brown, Member, Hydromechanics Committee, Hydraulics Divi- 
sion. 


Vibration Hydraulic Structures 


Campbell, Chief, Hydraulic Analysis Branch the Waterways 
Experiment Station. 


The Damping Oscillatory Waves Laminar Boundary Layers 


Eagleson, Assistant Professor Hydraulic Engineering, 
Massachusetts Institute Technology. 


Secondary Flow and Free Surface Phenomena Straight Open 
Channels 


J.W. Delleur, Associate Professor Hydraulic Engineering, 
Purdue University. 


Session Flood Control: 


Presiding: Eugene Fortson, Member, Executive Committee, and 


Williams, Jr., representing Flood Control Committee, Hydraulics 
Division. 

Problems Involved the Distribution Flow and Below the 
Latitude Old River 


Fred Chatry, Assistant Chief, Hydraulic Design Section, Corps 
Engineers, New Orleans, Louisiana. 


Sedimentation Aspects the Diversion Old River 


Fred Toffaleti, Chief, Hydraulic Model and Sedimentation Sec- 
tion, Mississippi River Commission, Vicksburg, Mississippi. 


Diversion from the Mississippi River Bayou laFourche, Louisi- 
and, Pumping Station 


Myers, Chief Engineer, Louisiana Department Public 
Works, Baton Rouge, Louisiana. 
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Maximum Results from Limited Available Information Hydrau- 
lic Computations 


Oakes, Chief, Hydraulic Section, Louisiana Department 
Public Works, Baton Rouge, Louisiana. 


Other Meetings 


Plans for the Reno Convention June 20-24, 1960, are progressing rapidly. 
The presiding Executive Committee member will Mr. Dickinson, Vice 
Chairman. The Sessions Tidal Hydraulics, Hydraulic Structures, and Flood 
Control will arranged Messrs. Ingerson, Harold Pratt, and 
Francis Christian, respectively. 

The papers and authors have not been confirmed, but the themes tidal 
hydraulics and hydraulic structures will Regime Influences Water 
Quality” and “Hydraulics High Head Outlet Works.” The theme the Flood 
Control Session not available this writing. 

Plans are being made for joint session with the United States Committee 
Large Dams the hydraulics large dams. More this the next 
Newsletter. 

Plans for the Hydraulics Division Seattle, Washington, are 
getting underway. Sessions are being planned Tidal Hydraulics, Hydrology, 
Hydromechanics, Sedimentation, Hydraulic Structures, and Flood Control. 

Preliminary plans for the Boston Convention, October 10-14, 1960, include 
sessions Tidal Hydraulics, Hydrology, Hydromechanics, Sedimentation, and 
Hydraulic Structures. 


Colorado Hydraulics Division Active 


The Hydraulics and Irrigation and Drainage Divisions the Colorado Sec- 
tion combine their meetings, each division being responsible alternate 
months for the meeting program. Some very interesting programs have been 
arranged both Divisions. One their recent and most interesting pro- 
grams included talk “The Blue Nile Water Resources Development 
Ethiopia,” Mr. Herbert Riesbol, Chief, Hydrology Branch, Division 
Project Investigations, Bureau Reclamation, Their future 
plans include talk Use Soil Moisture Relationships Planning 
Irrigation System,” Mr. Orville Parsons, Soil Scientist, Soil Conserva- 
tion Service, Lamar, Colorado. 

The officers the Hydraulics Division the Colorado Section for 1960 are: 


Chairman: Howard Matthai, Geological Survey, Building 25, Denver 
Federal Center, Denver, Colorado. 


Vice Chairman: Claude Hunter, Soil Conservation Service, Broadway 
and Belleview, Littleton, Colorado. 


Secretary-Treasurer: William Simmons, Jr., Bureau Reclamation, 
Building 56, Denver Federal Center, Denver, Colorado. 


Interesting Program 


Many ASCE members will interested the Hydraulics Division ASME 
Conference held Houston March 7-9, 1960. The program will include 
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two special symposia, one cavitation highlighting the problems face and 
research needs rather than specific accomplishments. The other stall 
and unsteady flow fluid machinery. 


For Your Calendar 


ASCE Meetings 


March 7-11, 1960 ASCE, New Orleans Convention 
June 20-24, 1960 ASCE, Reno Convention 
August 17-19, 1960 Hydraulics Division Conference, Seattle 
October 10-14, 1960 ASCE, Boston Convention 
April 10-14, 1961 ASCE, Phoenix Convention 
October 16-20, 1961 ASCE, New York Convention 
February 1962 ASCE, Houston Convention 
May 1962 ASCE, Omaha Convention 


October 15-19, 1962 ASCE, Detroit Convention 


Non-ASCE Meetings 


February 12-13, 1960 Regional Meeting, Society American 
Military Engineers, 
Jacksonville, Florida 

March 7-9, 1960 Gas-Turbine Power and Hydraulic Conference, 


American Society Mechanical Engineers, 
Houston, Texas 


April 4-8, 1960 Nuclear Congress and Exhibit, 
Coliseum, New York, 
April 28-30, 1960 Annual Meeting, American Geophysical 
June 18-21, 1962 Fourth National Congress Applied Mechanics, 


University California, Berkeley, California 


Use the Hydraulics Division Newsletter 


You are urged continue use the Division Newsletter for announce- 
ments, inquiries, personnel news, committee reports, surveys and other items 
interest Division members. short note summarizing the highlights 
committee meetings particularly requested. Suggestions for improvement 
the Newsletter will appreciated. 


Deadline dates for Newsletter contributions: April 1960 issue--February 20; 
June 1960 issue--April 20. 


James Ball, Editor 
Box 7416 
Lakewood 15, Colorado 
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